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ABSTRACT 


/ 

! 

\ 

> 

Certain  measurements  of  the  electrical  and  mechanical 
properties  of  transducers  are  of  importance  to  designers. 
These  can  be  deterained  froa  accurate  coaplez  admittance  and 
impedance  aeasnreaent3  made  throughout  the  frequency  range 
of  interest.  Sanaa  1  collection  of  needed  data  is  a 

tiae-consuaing  process  that  Is  prone  to  error.  The 
coaputerized  system  herein  described  substantially  reduces 
the  tiae  reguireaent  and  produces  more  accurate  output  than 
can  be  obtained  utilizing  manual  aethods.  In  addition  to  a 
general  discussion  of  equivalent  circuits  and  the 
instruments  employed,  this  report  includas  samples  of  plots 
and  calculated  parameters  for  piezoelectric  and 
aagnetcstrictive  transducers,  the  experimental  comparison 
with  traditional  manual  aethods  and  the  program  listing  for 
two  versions  of  the  system  allowing  varying  amounts  of 
operator  options. 
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Transduction  is  defined  ia  an  extremely  general  sense  as 
any  conversion  of  one  form  of  energy  into  another.  The 
devices  responsible  for  this  conversion,  in  either 
direction,  are  transducers.  For  this  author,  the  interest 
is  the  conversion  of  alectricity  into  sound  and  vice  versa. 
Comeonly,  a  transducer  that  converts  electrical  energy  into 
sound  is  called  a  projector,  loudspeaker,  source  or 
transmitter .  Those  serving  the  reverse  function  are 
coeaonly  called  hydrophones,  licrophones,  or  receivers. 
Some  types  of  transducers  are  reversible  and  are  capable  of 
filling  either  role.  This  study  is  restricted  to  the 
analysis  of  reversible  transducers. 

&.  SECHAHISas  OF  TfiiH  SDU3TI0N 

There  are  two  basic  types  oE  mechanises  commonly  used  in 
reversible  electroacoustic  transducers  designed  for 
underwater  sound,  one  involving  interaction  between 
mechanical  motions  and  forces  and  electrostatic  fields  and 
the  other  between  seohanical  lotions  and  forces  and  magnetic 
fields.  Both  types  are  aaeaable  to  analysis  using 
equivalent  electrical  circuits;  but,  due  to  the  basic 


differences  in  conventions  used  to  describe  forces  in 
electric  and  in  magnet ic  fields,  the  details  are  different. 
Magnetically  coupled  transducers  of  the  soring  coil  type 
have  the  diaphraa  attached  to  a  wire  coil  suspended  in  a 
Magnetic  field.  It  i3  the  interaction  of  the  current  in  the 
coil  and  the  magnetic  field  that  induces  a  force  on  the  coil 
and  causes  the  diaphraa  to  aove  and  radiate  an  acoustic 
wave  [Ref.  1].  These  classifications  say  be  further 
subdivided  based  on  the  sotor  mechanises  involved  (i.e. 
electrodynamic,  electrostatic,  aagnetic,  aagnetostrict ive, 
and  piezoelectric)  [Ref.  2].  The  classification  systems  for 
transducers  vary  depending  upon  the  analysis  being  done. 
Descriptions  of  various  classifications  are  presented  by  the 
previous  references  as  well  as  by  Heaslip  [Ref.  3],  and 
Sherman  [Ref.  4  ]. 

B.  PURPOSE  OF  THIS  3T0DY 

411  electroacoustic  transducers  may  be  considered  as 
being  composed  of  eLectrical  elements  so  that  the  analysis 
of  the  transducer  is  greatly  simplified  if  an  equivalent 
electrical  circuit  is  used. 

While  the  most  common  transducers  currently  in  use  in 
underwater  sound  are  made  of  piezoelectric  materials,  the 
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aim  of  this  study  is  to  provide  a  generalized, 
semiautomat  ed  means  of  obtaining  aeaningful  transducer 
aeasureaents  on  any  type  of  transducer  with  the  operator 
having  mini  oil  knowledge  of  the  device. 

The  equipaent  set-up  utilizes  basic  devices  (e.g. 
voltmeter,  sweep  oscillator,  Dranetz  impedance  meter, 
aaplifier)  interfaced  to,  and  controlled  by,  a  desk-top 
computer.  This  interfacing  provides  eore  precise  data 
acquisition  than  can  be  obtained  from  manual  plots 
interpreted  by  an  operator.  The  time  saving  factor  is  a 
major  motivation  for  this  choice. 

Certain  measurements  of  the  electrical  and  mechanical 
properties  of  the  transducer  are  of  importance  to  designers. 
In  some  cases,  it  is  desireable  to  have  the  resonant 
frequency  fall  within  a  particular  range.  More  likely, 
there  will  be  concern  for  a  high  electromechanical  coupling 
coefficient.  Efficiency  may  be  the  major  concern  or  perhaps 
it  is  desired  to  have  high  losses  and,  therefore,  a  low 
mechanical  quality  factor  with  consequent  "flatness"  of 
response  [Bef.  5j. 

Direct  measurement  of  these  values  is  not  always 
possible.  However,  many  of  the  electrical  component  values 
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needed  to  calculate  lechanical  parameters  can  be  readily 
obtained  frca  propacLy  annotated  iapadance  and  adaittance 
plots  for  a  particular  transducer. 

Currently,  theca  are  automated  systems  designed  to 
analyze /transducers,  such  as  the  wqh-12  [Ref.  6].  Pew  are 
saall  enough  to  be  readily  transported  to  the  transducer  as 
is  this  systee. 

C.  POR1UT  OP  THE  REPOST 

The  basic  physical  and  electrical  theory  used  in 
preparation  of  this  report  is  presented  in  chapter  II. 
While  aore  thorough  analyses  of  transducers  and  equivalent 
circuits  may  be  found  in  some  of  the  references,  only  data 
pertinent  to  the  measurements  obtained  by  this  system  are 
presented. 

Chapter  III  erpLains  the  interfacing  between  various 
components  of  the  equipment  set-up  and  the  computer.  The 
operator/computer  interaction  is  also  addressed. 
Representative  samples  of  computer  output  are  included. 

Chapter  17  briefly  describes  the  testing  and  evaluation 
of  the  system  and  program.  The  results  and  conclusions  are 


presented  here 


The  cob pater  progress  erittaa  in  BASIC  computer  language 
are  included  in  Appendix  A. 
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Sines  such  of  the  physics  and  engineering  involved  in 
discussions  of  transducers  is  common  to  both 
electric-coupled  and  magnetic-coupled  transducers,  this 
report  will  sake  a  general  approach  and  point  out 
differences  where  thay  occur.  Table  I,  at  the  end  of  this 
chapter,  contains  a  listing  of  symbols  used  in  this  report 
and  refers  to  symbology  usad  in  some  of  the  references 
cited. 

1.  GENERAL 

The  approach  usually  taken  in  electroacoustic  analysis 
is  to  use  an  electric  circuit  :d  model  the  transducer.  In 
simplest  fora  one  may  use  a  two-port  network  with  one  port 
representing  the  electrical  terminals  and  the  other  the 
mechanical  terminals.  Analogs  for  force  and  speed  may  be 
voltage  and  current  or  vice  versa  depending  on  the  type  of 
transducer. 

The  canonical  equations  describing  the  behavior  of  the 
two  pore  transducer  may  be  written  as: 

7»Z0-  I  ♦  Tem-0  (2.1) 

and 
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(2.2) 


♦  z«-° 

where  7  and  I  are  the  voltage  and  current  at  the  electrical 
port  and  F  and  0  are  the  force  and  velocity  at  the 
aechanical  port.  Iem  is  a  transduction  coefficient  relating 
the  electromotive  force  at  the  electrical  port  to  the 
velocity  in  the  aechanical  net,  and  Tme  is  a  transduction 
coefficient  that  relates  the  force  developed  at  the 
aechanical  port  of  a  two-port  network  to  the  current  in  the 
electrical  aesh  [Bef.  8],  <i>  is  the  transformation  factor 

and  aay  be  defined  as  <b=T*m/Z n.  The  canonical  equations  for 
a  maqnetically  coupLed  transducer  of  the  moving  coil  type 
are: 


F-  -B •  1  •  I  ♦  Zm3  (2.3) 

7-Z0-I  ♦  B •  1  •  0 

where  3  is  the  magnetic  field  and  1  is  the  length  of  the 
wire  in  the  coil.  (Mote:  This  choice  is  valid  for 
eiectroayna aic  transducers  (aoving  coil)  as  written  using 
B'1*Ten»*  Magnetostcictive  transducers  obey  the  same  set  of 
equations  but  B-l.) 
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For  el ectromechanical  coupling.  Figure  2.1  shows  typical 


circuits  and  transf orsations  'Bef.  7]. 


iraiBi; 


Pig.  2.1.  Typical  Circuits  and  Transformations 

In  the  vicinity  of  resonance,  which  is  of  concern  here, 
the  Mason  circuit  or  aodified  Via  Dyke  circuit  is  valid.  It 
is  also  valid  for  low  frequency  values  which  will  also  be 
needed  for  some  calculations  ‘  Ref.  8].  This  circuit  is 
shown  in  Figure  2.2  where  Lm  ,  Rm  ,  and  Cm  represent 
mechanical  values.  C0  is  the  clamped  or  blocked  or  shunt 
capacitance . 
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Figure  2.2.  (a).  Boiified  Van  Dyke  or  Bason  circuit  (b)  . 

Basic  Transducer  Circuit. 

Kinsler ,  and  others  [Bef.  1].  define  the  following 
measurable  mechanical  and  electrical  iapedances  of  a  system: 

Blocked  electrical  impedance  (ohms)  *Z0*V/IIu,1iq  _  51 


Free  electrical  impedance  (ohms) *2  “V/Il^g 


(2.6) 


Open  circuit  mechanical  iapsdance  (N  •  s/m)  *Znf,F/lI|j3Q 

(2.7) 


Short  circuit  mechanical  Lmpedance  (8-  s/a)  ^Z^^P/g  iy_Q 

(2.8) 

Magnetically  coupled  transducers  say  be  represented  by 
the  circuit  of  Figure  2.3,  where  4>m  is  a  transform tion 
factor.  4»m  *  B*1  for  moving  coil  transducers  [8ef.  2]. 


Figure  2.3.  Equivalent  Circuit  for  a  Magnetically  Coupled 
Transducer. 

Important  properties  of  the  transducer  may  be  found  by 
studying  the  electrical  impedance  of  the  system  as  a 
function  of  driving  frequency.  The  unloaded  driving  point 
electrical  impedance  may  be  determined  from  the  canonical 
equations  (2. 1-2.4)  aid  is  given  by 

z-vi^o-z°  ♦  <-Wme/zmi  <2-9> 
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Similarly,  the  short  circuit  mechanical  driving  point 
impedance  is  given  by 

Zmijj* F /U )y_Q  *Zm  ♦  (2.10) 

The  electrical  impedance  is  influenced  by  the  motion  of  the 
coupled  mechanical  system,  as  indicated  by  the  second  term 
on  the  right  hand  side.  This  term  is  referred  to  as  the 
motional  impedance,  and  defined  by  Hunt  [Bef.  2],  as 

zmota  (-WWH/y  (2.11) 

Similarly,  the  mechanical  driving  point  impedance  is 
influenced  by  the  electrical  load. 

B.  USE  OF  COMPLEX  IMP  ED AN CE/A OMITTANCE  PLOTS 

1-  IllS  Resosaase  SlKla 

Graphical  displays  of  iapedance  aid  admittance  data 
can  be  extremely  useful  in  the  analysis  of  a  transducer. 
When  properly  scaled  and  annotated,  quantitative  results  may 
be  obtained  from  the  diagrams.  generally,  in  the  vicinity 
of  resonance,  the  element  is  regarded  as  having  a  single 
degree  of  freedom  and  the  locus  of  points  of  the  complex 
iapedance  is  a  circle  in  the  neighborhood  of  the  resonance 
frequency  [Bmf.  7],  [Ref.  3],  Figure  2.4  and  Figure  2.5 
are  typical  resonance  circle  diagrams  where,  in  Figure  2.4, 
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Fig.  2.5.  Complex  Maittance  Diagram 

In  Figure  2.5,  f2  and  f3  equate  to  the  frequencies 
of  the  half  power  points.  F 0  is  the  frequency  of  resonance. 
The  dashed  line  would  be  the  same  ideal  transducer  tested 
under  load  conditions.  (Note:  Susceptance  may  be  entirely 
negative  for  a  aagnetically  coupled  transducer.) 

Motional  impedance  components  may  be  obtained  using 
the  following  formulae: 
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®mot  *®E  (2.12) 

*mot  **E  +  3^0  (2.13) 

with  Zjjjqj  *  8mo^  *  j  ^mot  (2.14) 

The  motional  admittance  components  are  given  by: 

Gmot  3  ge  •  <1/Ro>  (2.15) 

Bmot  3  ®e  ”  ^  J  C0  (2  •  16) 

with  Tmot  *  3mot  ~  j  ®mot  (2.17) 


The  frequency  of  maximum  power  output  at  constant  voltage 
input  is  called  the  aechaaical  resonance  and  is  represented 
in  Figure  2.6  by  f9  [Ref-  9]. 

In  the  motional  admittance  diagram.  Figure  2.7,  fO 
is  the  frequency  of  maximum  conductance  and  resonance.  F2 
and  f3  are  the  half  power  points. 

2.  £2SL20fl®£!  SiEiOHIii 

Impedance  measurements  are  made  at  the  electrical 
terminals  of  the  transducer.  The  behavior  of  the  real  and 
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Fig.  2.6.  Notional 

Impedance 

Diagram 

imaginary  components  of  the 

complex 

impedance 

may  be 

observed  as  the  frequency 

spectrum 

is  swept 

at  a 

sufficiently  slow  rate  so 

that 

the  syst 

em  is 

quasi-steady-state.  Hhen  the 

n  echanical 

system  is 

cla  mped 

so  that  it  cannot  move,  the  blocked  electrical  impedance  nay 
be  measured  and  used  to  calculate  notional  impedance,  since 
it  is  difficult  to  adequately  claap  a  device  over  the  entire 
frequency  spectrum,  values  nay  be  obtained  far  above  and 
below  resonance  and  the  intermediate  values  inferred.  This 
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Fig.  2.7.  notional  hdmittance  Diagram 

is  especially  important  for  tae  reactance  and  susceptaace 
curves  whose  clamped  values  vary  more  drastically  with 
frequency  than  do  the  real  components  of  admittance  or 
impedance. 

The  maximum  walue  of  motional  impedance  will  equate 
to  the  mechanical  resistance  and  is  the  diameter  of  the 
circle.  This  value  will  be  of  importance  in  calculating  the 
efficiency  of  the  system. 
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Plots  of  various  input  electrical  parameters 
(conductance  (GE)  and  susceptarca  (Bg  )  or  resistance  (H£) 
and  reactance  (I  E )  versus  frequency)  can  provide  all 
of  the  information  necessary  to  obtain  the  electrical  and 
■echanical  quality  factors  for  the  systea.  These  values 
indicate  the  sharpness  of  the  resonance, 
a.  Impedance 

The  circuit  of  Figure  2.8  is  useful  for  analysis 
of  the  behavior  in  the  vicinity  of  resonance. 


Pig.  2.8.  Equivalent  Circuit  Diagraa 

A  plot  of  the  input  electrical  impedance  components  versus 
frequency  will  resemble  Figure  2.9. 
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Co  B  VO 


Fig. 


2.9.  Besistanca  and  Beactaace  Plotted  Versus  Frequency 


b.  idnittanra 

The  aotioaal  adaittance  equations  are 


W  *  smot  *  jB 


mot 


_  1*mot-jXmot_ 


®niot+j^mot  ®mot 


®mot 


—  ^mot 


®mot  *  ^mot 


®mot  *  ^mot2 


(2.18) 


Plots  of  conductance  and  susoeptance  versus  frequency  for 
Figure  2.8  will  reseible  Figure  2.10.  The  frequencies  where 
the  susceptance  (8)  is  zero  correspond  to  the  electrical 
resonance  and  antiresonance  frequencies. 


Figure  2.10.  Conductance  and  Susceptance  Plotted  Versus 
Frequency 

C.  EVALUAT ION  OF  T  RAN  SDOCER  PARAMETERS 

This  systea  is  designed  to  obtain  useful  inforaation  on 
devices  used  in  transacting  ani  receiving  in  water.  As 


such#  it  is  first  necessary  to  obtain  data  in  an  unloaded 
situation  in  order  to  be  able  to  separate  the  effects  of 
acoustic  loading.  for  transducers  designed  for  use  in 
water,  air  noraally  provides  an  adequate  medium  for 
evaluation  cf  the  acoustically  unloaded  properties.  Placing 
the  transducer  in  water  aeans  the  two-port  network  will  now 
reflect  the  sum  of  the  aechaaical  and  the  load  impedances. 
Then,  according  to  Hunt  [fief.  2], 

ZT  -  Zm*  2,  *  <2*19> 

or 

Zm  *  z,  *  <Rm  *■  R  i )  ( 1  ♦  j  ‘  2  •  Q  •  p)  (2.20) 

where  p  is  a  frequency  parameter-- (p®. 5  ( (^4j0  )*  (u;0Aj)  )  and  a 
is  the  quality  factor  measured  for  the  resonance  (i.e. 
electrical  or  mechanical  depending  on  the  type  of 
measurements) . 

Data  obtainable  from  admittance  asasurements  will 
provide  needed  values  for  determining  the  motional 
admittance  circle.  The  maximum  value  of  conductance  will 
occur  at  the  frequency  of  electical  resonance.  The  quality 
factor  can  be  determined  from  the  frequencies  at  which  the 
conductance  assumes  half  the  iiffference  of  its  maximum 


value  and  its  blocked  value  (i.  e.  Q  =  f  /(fu  -fj).  These 
frequencies  sill  also  coincide  with  the  frequencies  (f2  and 
£3)  of  the  local  maximum  and  minimum  susceptance  values  near 
resonance.  Blocked  capacitance  may  be  calculated  from 
measurements  of  susceptance  far  above  and/or  below 
resonance.  For  an  electrically  coupled  transducer  with 
small  dielectric  losses,  valid  measurements  of  CQ  may  be 
made  very  far  b9low  resonanca.  For  magnetic  coupling, 
measurements  must  be  taken  far  above  and  far  below  resonance 
and  the  values  between  interpolated  linearly  between  these 
two. 

Impedance  data  will  provide  the  frequency  of 
antiresonance.  The  frequencies  of  the  half  values  for  the 
difference  between  the  maximum  resistance  and  the  blocked 
resistance  will  correspond  to  tae  maximum  and  minimum  values 
of  reactance  and  allow  for  calculation  of  a  mechanical 
quality  factor.  Blocked  capacitance  or  inductance 
(depending  on  the  type  of  transducer)  may  be  calculated  from 
reactance  data  obtained  far  above  and/or  below  resonance. 

The  following  paragraphs  briefly  describe  the  definition 
and  calculation  of  various  constants. 
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1.  fiuaOi££  &2i2E  z  2 


This  factor  indicates  the  sharpness  of  resonance. 


frequency  of  aaz.  3  or  R 

freq.  difference  of  min. 
and  sax.  B  or  X 


f  (resonance) 
f  (upper)  -f  (lower) 

(2.27) 


2.  sisst£2assh§aical  ssaeiisa  saslfisisai.  =  ii(eSf) 

This  coefficient  indicates  how  tightly  the 
electrical  and  mechanical  meshes  are  coupled  [Ref.  2]. 

K2  (eff)  *  1-  (f B/f 0)  2  for  electromagnetic  coupling 

(2.22) 

and 

K*(eff)  *  1-(f3/f8)2  for  electrostatic  coupling 

(2 .23) 

where  f8  is  used  to  indicate  the  frequency  of  maximum 
resistance  and  fD  indicates  the  frequency  of  maximum 
conductance. 

3-  stasis  SaislUa.  £2&££UU&1  z  SI 

This  term  is  the  ratio  of  the  stored  mechanical 
energy  to  the  total  energy  stored  in  the  device  [Ref.  1  ], 
[Ref.  5]. 
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(2.  24) 


Co  +  Cj 


4.  at aat  £&£&£Uiasa  *  £o 

This  is  the  blocked  capacitance  of  the  system. 
There  is  nc  simple  *athod  of  measuring  this  quantity  (with 
the  possible  excaption  of  instruments  designed  far  a 
specific  frequency) .  Heasucaaents  must  be  made  at  a 


fr equency 

far 

re  mov  ed 

from 

frequencies 

of  mechanical 

resonance 

in 

order 

to 

minimize  the 

effects  of 

electromechanical  coupling.  If  the  capacitance  ratio 
(Cj/C0)  is  greater  than  fifty,  the  inaccuracy  should  be 
minimal  and  may  be  cansidarad  negligible. 

1  3 

CD  * - -  — —  where  &>  <<  <u0  (z  .  25) 

OjXo  Cl* 

5.  £liJIESd  Sssiilansi  *  tfl  a£.d  la 

R o  and  X0  may  be  closaLy  approximated  by  obtaining 
data  far  below  and  above  resonance  and  interpolating  between 
since  there  is  a  positive  frequency  dependence.  However, 
for  many  practical  purposes,  (and  for  electrically  coupled 
transducers),  data  may  be  obtained  far  below  resonance  which 


34 


provides  a  reasonable  approximation  to  the  actual  value 
throughout  the  frequency  range  of  interest.  These  values 
are  assuaed  to  represent  those  which  would  be  measured  if 
the  aechanical  parts  of  the  system  were  blocked  or  clasped 
to  prevent  notion. 

6.  aiaslsfi  x&diu&ii&s.  s  ko 

This  is  the  inductance  of  the  blocked  system 
(usually  with  magnetic  coupling!  and  is  obtainable  from  data 
collected  far  below  and  above  resonance. 


L 


X(interpola  ted) 


(2.26) 


7.  Cap acit ancs  Siiio  = 

This  is  the  ratio  of  the  blocked  capacitance  (C0)  to 
the  motional  capacitance  (C,). 

PI  (1-K2)  where  K2  is  the 

- -  static  coupling 

8  S 2  coefficient 

(2.27) 
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8.  EUitti  fi£  M§iit 


The 

figure  of  merit 

is  little 

mentioned  in 

the 

references 

probably  because 

it  is  an 

arbitrary  formula 

depending 

upon  a  designer's 

criteria. 

Miller  [Ref. 

5]r 

mentions  K2 

•Qm  as  one  possibility  and 

utilization  of 

the 

coupling  coefficient  as  another.  No  specific  calculation  of 
this  value  is  done  ia  this  report,  leaving  it  to  the  user  to 
insert  as  desired. 

9.  fliaaafcaiS 

Impedance  or  admittance  circles  ace  transformed  to 
motional  diagrams  by  subtracting  blocked  values  or  blocked 
values  multiplied  by  angular  frequency  as  discussed  in 
section  fl. 1  above.  I  he  diameters  of  the  resulting  circles 
will  provide  needed  values  for  efficiency  calculations.  D^, 
and  signify  the  diameters  of  the  motional  circle  for  the 
transducer  in  water  and  in  air,  respectively. 

D.  EPFICIE  NCI 

The  overall  efficiency  of  a  transducer  is  the  ratio  of 
the  power  delivered  to  an  external  load  connected  at  the 
output  terminals  to  the  total  power  at  the  input. 

At  resonance,  this  can  be  easily  calculated  from  known 
values.  For  an  electrically  coupled  transducer,  the 
efficiency  is  given  by: 
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where  and  T)w  are  the  diameters  of  the  motional  admittance 
circles  and  where  3  is  the  conductance  at  resonance. 

Similarly,  the  efficiency  at  resonance  for  a 

magnetically  coupled  transducer  is  given  by: 

D,.,(Dj  -D„.) 

E  *  1 (2.29) 

R 

where  these  are  diameters  of  the  motional  impedance  circles 
and  R  is  the  resistance  at  resonance  [Hef.  10]. 

The  mechanical  power  utilization  factor  is  part  of  the 
overall  efficiency.  It  is  given  by  Rj/fR^R^)  which  can  be 
expressed  in  terms  of  loaded  and  unloaded  values  of  Q  or  by 
the  diameters  of  circles. 


Rm  ♦  R/ 


(2.30) 


(Note:  All  values  should  be  taken  from  either  admittance 


data  or  impedance  data  but  should  not  be  mixed.) 


Once  the  efficiency  at  resonance  is  determined, 
utilization  of  the  impedance  data  of  a  magnetically  coupled 
transducer  allows  for  calculation  of  the  lagging  phase  angle 

<p>. 


P*  0.5-  ArcCos  ( (Rres  "Ro> /Dw)  (2.31) 

The  frequency  parameter  for  maximum  efficiency  may  be 
calculated  as 

p*  (DwSin(2-p)  > /(4-E0. 3wi  (2.32) 

and  from  the  earlier  relationship  between  p  and  the 
frequency  ratios,  tha  frequency  for  maximum  efficiency  f or  a 
magnetically  coupled  transducer  may  be  calculated. 

f  *  f0  (2'P  ♦  (p2  +  i)  i/fe)  (2.  33) 

The  frequency  of  maximum  efficiency  may  not  coincide 
with  either  the  resonant  frequency  or  the  frequency  cf 
maximum  admittance.  For  the  aagne+osfrictive  transducer, 
this  frequency  of  maximum  efficiency  is  usually  greater  than 
that  at  resonance. 
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A  development  by  Camp  *  Bef .  8]r  shows  that  a 
magnetostr  ictive  transducer  may  produce  more  acoustic  power 
at  the  resonant  frequency  even  though  efficiency  may  be 
less.  Normally,  the  frequency  for  optimum  acoustic  power 
will  fall  between  the  mechanical  resonance  frequency  and  the 
frequency  of  optimum  efficiency.  This  is  due  to  eddy 
current  losses  and  magnetic  hysteresis. 

Since  transduction  losses  are  low,  maximum  power  output 
for  piezoelectric  devices  is  limited  mainly  by  voltage 
limits.  The  optiium  operating  frequency  is  accepted  as  the 
frequency  of  maximum  admittance  *  Ref .  8]. 

If  one  assumes  that  the  transducer  has  a  high  Q  (sharp 
resonance)  and  that  the  blocked  resistance  is  constant 
throughout  the  frequency  range  for  the  impedance  circle, 
then  the  potential  efficiency  may  be  expressed  as 


Pot.  Eff. 


(2.34) 


The  potential  efficiency  represents  the  maximum  efficiency 
available  for  the  most  favorable  loading  conditions.  This 
equation  deletes  the  requirement  for  data  from  a  loaded 
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condition  but  requires  the  air  motional  impedance  loop  to 


very  closely  approximate  a  circle  [Ref.  9]»  "High  Q w  as 
used  here  implies  values  greater  than  thirty.  If  the  air 
■otional  impedance  loop  is  not  approximately  circular  and/or 
the  Q  is  loe,  then,  from  reference  2, 

pot.  Eft.  -  J!.J 4s2Uii5-:-ii|.!SEAlS~ 

(  1  {£  COS 2  0  >  V*  ♦  (1-  jg  SIN*  p  ) 1/2 

{2.35) 


E.  SUMMARY 

There  is  helpfuL  information  to  be  obtained  from  both 
impedance  and  admittance  data.  Some  transducers  more 
readily  reveal  properties  through  one  type  of  diagram  than 
another,  Normally,  admittance  data  is  collected  for 
electrically  coupled  systems,  voile  impedance  data  is  better 
for  studying  magnetically  coupled  systems  [Ref.  2]%  Due  to 
the  ease  of  obtaining  data  with  the  system  proposed  here, 
both  sets  of  data  and  diagrams  say  be  rapidly  obtained. 


'4.0 


TABLE  I 

COMPARATIVE  SYMBOLOGY 
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ckt- 
Mech.  Imped. 

Mechanical 

Resistance 

Mechanical 

Capacitance 

Motional 

Impedance 

Motional 

Admittance 

Motional 

Resistance 

Motional 

Reactance 

Motional 

Conductance 

Motional 

Susceptance 

Mechanical  , 

Resonance 

Electrical  ( 

Resonance 

Mechanical 
Quality  Factor 

Electrical 
Quality  Factor 
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ill.  A24£XMI2fi  12  l  :2H£2IMIZ52  mna 

In-phase  and  out-of-phase  components  of  the  admittance 
or  impedance  are  needed  in  the  analysis  of  a  transducer.  An 
initial  attempt  was  lade  to  usa  a  lock-in  analyzer,  but 
ultimately  a  Dranetz  impedance  aeter  was  used. 

A.  EQUIPMENT  SET-OP 

Data  needed  foe  the  calculation  of  properties  of 

transducers  may  be  obtained  using  a  lock-in  analyzer.  The 
magnitude  and  phase  or  in-phase  and  quadrature  components 
will  indicate  resonances  and  provide  the  necessary 
information. 

The  initial  set-up  utilized  the  Princeton  Applied 

Research  Company  Model  5204  Lock-in  Analyzer.  Phase 

information  would  have  allowed  for  desired  calculations.  In 
experimentation,  the  circuit  required  for  obtaining  valid 
phase  information  became  too  coaplex  to  be  practical.  The 
lock-in  analyzer  is  designed  to  operate  in  the  1  Hz  to  100 
KHz  frequency  range  although  it  was  determined  to  be 

accurate  at  higher  frequencies.  At  low  frequencies,  there 
is  a  low  source  resistance  ar.d  the  signal  to  noise  ratio 
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becaae  too  small  to  provide  &  reference.  It  was  concluded 
that  coupler  iapedaare  measures ents  would  be  another  means 
of  getting  the  necessary  data. 

As  briefly  mentioned  in  the  introduction,  there  are  many 
methods  for  obtaining  impedance  information.  Since  this 
au+hor  desired  a  portable  system  usable  for  underwater 
transducers,  a  bridge  network  seemed  reasonable. 

The  set-up  used  to  design  this  system  is  shown  in  Figure 
3.1.  The  Hewlett-Packard  HP-35  computer  directs  the 
peripherals  for  the  data  taking.  The  printer  and  plotter 
are  non-essential  to  the  system  and  were  used  to  provide 
larger  format  graphics  and  printout  than  are  available  with 
the  built  in  thermal  printer  provided  with  the  HP-85. 

1.  Hewlett -Packer d  H£-§.5  Computer 

The  HP-85  is  an  eight  bit  microcomputer  that 
utilizes  BASIC  computer  language.  It  has  a  127  millimeter 
diagonal  black  and  white  electromagnetic-deflection  CRT.  A 
32  character  per  line  thermal  printer/plotter  is  part  of  the 
unit.  The  computer  has  16K  bytes  of  read/write  memory  of 
which  14,  579  are  available  to  the  operator.  For  this 
program  a  memory  module  is  needed  to  expand  the  memory  to 
32K  bytes.  Programs  or  data  may  be  stored  on,  or  read  from. 


magnetic  tape  cartridges.  To  interface  with  peripheral 
equipment,  an  I/O  ROM  and  an  interface  card  were  added  to 
provide  HP-IB  (IESE-438)  instrumentation  interface 
capabilities. 

2.  §Y4thSsj.Z3C/?1» action  isiar^tor 

The  Hewlett-Packard  Model  3325A  Synthesizer/Function 
Generator  can  produce  five  different  continuous  waveforms. 
For  this  application  a  sine  wave  was  desired  for  which  this 
instrument  has  a  frequency  range  of  1  microhertz  to  20 
megahertz.  Frequency  may  be  specified  with  up  to  eleven 
digits  of  resolution.  Output  amplitude  is  7  millivolt  to  10 
volts  peak-to-peak  into  a  fifty  ohm  load.  This  nodal  is 
fully  programmable  through  the  rear  panel  Hewlett-Packard 
Interface  Bus  (HP-IBJ  . 

3.  Dr a  net 3  Complex  Impedance — Admittance  Meter 

The  Model  103C  Complex  Impedance-- Admittance  Meter 
(CIAM-100C)  is  a  transistorized  instrument  that  measures 
vector  impedance,  vector  admittance,  vector  amplitude  and 
phase.  It  covers  a  frequency  range  from  100  Hz  to  20D  KHz 
in  three  ranges  of  3.  1  to  2  KHz,  1  to  20  KHz  and  1 0  to 
200RHz.  For  measuring  impedance  a  constant  current  is 
produced  in  the  unknown  impedance. 


The  amplitude  of  the 


Figure  3.1  Equipment  Set-up 


currant  is  set  by  the  range  switch.  The  voltage  developed 
across  the  unknown  iipedance  is  amplified  and  fed  to  two 
resolvers.  Por  measuring  admittance,  a  constant  voltage  is 
applied  to  the  unknown.  The  current  through  the  unknown  is 
measured  by  means  of  voltage  developed  across  a  small 
range- selected  resistor  placed  in  series  with  the  unknown. 
The  resolver  outputs  are  proportional  to  conductance  and 
susceptance  of  the  unknown.  The  claimed  accuracy  of  the 
instrument  is  plus  or  minus  two  percent  full  scale 
amplitude.  The  Orauetz  meter  requires  an  input  signal 
voltage  of  3  Vrms.  In  taking  measurements  on  an  unknown 
transducer,  initially  the  frequency  of  interest  may  not  be 
known  as  well  as  the  appropriate  scale  settings  for  the 
meter.  Measurements  made  on  the  Dranetz  meter  using  a 
variable  resistance  box  and  voltmeters  indicate  linearity 
throughout  the  range  settings.  However,  it  is  necessary  to 
insure  measurements  do  not  cause  either  meter  range  to  be 
exceeded.  The  limit  indicator  lights  are  activated  when  the 
unknown  impedance  or  admittance  is  substantially  greater 
than  the  full  scale  setting.  ihen  this  occurs,  an  internal 
relay  shunts  the  meters  and  the  DC  output  terminals  to  zero. 


4.  2a ta  2ai£ 

The  Hewlett-Packard  Model  3497A  Data 
Acquisition /Control  On  it  is  used  here  as  a  voltmeter  and 
scanner  to  aeasure  the  DC  output  voltage  proportional  to  the 
deflection  of  the  maters  of  the  Dranetz.  The  instrument 
measures  and  displays  voltaga  values  to  five  and  one  half 
digits.  It  is  internally  triggered  by  software  command 
during  the  data  taking  but  takes  the  next  reading  as  soon  as 
the  first  is  completed  for  obtaining  averages. 

5.  £o we£  £m£liSiiI 

This  instrument  is  not  assential  to  the  operation  of 
the  system.  The  Dranetz  metar  is  designed  to  accept  direct 
ac  voltage  input;  however,  the  particular  meter  used 
exhibited  more  stability  when  the  amplifier  was  included  in 
the  set-up.  The  amplifier  had  unity  gain  and  was  only  used 
to  lower  the  effective  output  impedance  of  the  frequency 
synthesizer  from  fifty  ohms  to  two  milliohms. 

6.  eiicussiaa 

This  equipment  set-up  is  suitable  for  an  initial 
system.  A  system  designed  for  regular  and  routine  use 
should  include  a  different  impedance  measuring  system  than 
the  Dranetz.  Because  of  the  design  feature  of  the  Dranetz, 
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a  great  deal  of  operator  interface  is  required  to  change  the 
various  different  scale  and  range  factors  which  cannot  be 
set  directly  by  the  coeputar.  It  is  too  easy  for  an 
operator  to  miss  changing  oaa  setting  or  to  overlap 
frequency  range  scales  in  the  spectrum  sweep  and,  thereby, 
obtain  erroneous  data.  Also,  the  scale  factors  aust  be 
manually  entered  prior  to  the  computer  doing  any  cal¬ 
culations.  The  Dranetz  also  had  annoying  features,  such  as 
drift  of  the  zero,  excessive  noise  on  low  impedance  ranges, 
and  DC  output  at  full  scale  deflection  of  0.8  -  0.95  volts. 

B.  DISCUSSION  OP  THE  PROS  BAM 

The  HP-85  is  the  controlling  unit  directing  the  other 
instruments  to  take  or  send  dati  as  needed.  Memory  storage 
is  handled  by  this  mic rocoeputer . 

The  goal  was  to  design  a  program  than  would  minimize  the 
time  necessary  to  collect  pertinent  data  while  maximizing 
accuracy.  The  program  was  to  accommodate  different  media, 
allow  for  collection  of  admittance  or  impedance  data  or 
both,  and  work  for  transducers  with  either  electric  or 
magnetic  coupling,  while  being  easy  to  operate.  This  was 
accomplished;  but  due  to  the  large  differences  in  properties 
of  different  transducers,  a  second  program  was  written  tc  be 
used  by  more  experienced  operators. 
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One  program  is  designed  so  that  a  person  with  minimal 
knowledge  of  the  transducer  or  measurement  procedures  may 
run  it  and  obtain  meaningful  information.  It  is  suggested 
that  the  operator  be  familiar  with  the  operation  of  the 
interfacing  instruments  (especially  the  Dranetz  meter).  An 
equivalent  program  is  designed  so  that  the  operator  may 
select  the  bandwidth  for  data  collection  and  may  specify  the 
frequency  fcr  measurement  of  tie  blocked  values  (done  as  a 
percentage  of  the  resonant  frequency) . 

The  prograa  is  broken  into  three  sections.  The  first 
gives  the  operator  an  overview  of  the  conductance  and 
resistance  of  the  transducer  uider  test.  Fart  two  is  the 
data  collection  for  desired  measurements  with  options  for 
plots  and  lists  of  the  data.  The  final  section  performs  the 
calculations  for  motional  data,  efficiencies,  and  other 
desired  outputs.  The  following  subdivisions  address  the 
program  in  detail.  Figure  3.  2  is  a  flow  chart  of  the 
program.  The  programs  written  ir.  BASIC  are  included  as 
Appendix  A. 

This  first  portion  has  been  included  to  allow  for 
visual  and  graphic  display  of  tae  response  of  the  transducer 
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over  a  selected  frequency  range.  The  operator  is  required 
to  insert  information  on  the  type  of  coupling  (magnetic  or 
electric)  of  the  transducer  and  the  medium  (air  or  water)  in 
which  measurements  are  to  be  taken.  The  operator  is 
instructed  to  set  the  Oranetz  for  the  measurement  of 
admittance  data  and  to  specify  the  range  to  be  cowered  in 
the  frequency  sweep.  The  opentor  also  enters  the  desired 
voltage  to  be  sent  to  the  Dtanetz.  This  manual  voltage 
input  is  necessary  since  at  higher  frequencies  lower  rms 
voltages  are  necessary  to  keep  the  a c  voltage  meter  on  the 
Dranetz  at  3.0  V. 

Through  the  Data  Acguisition  Unit,  the  relative 
amplitude  of  the  real  part  of  the  admittance  is  sampled  at 
300  equally  spaced  points  throughout  the  selected  frequency 
range.  k  graph  is  displayed  on  the  CRT  and  then  copied  on 
the  thermal  printer  of  the  HP-35  for  later  reference.  In  the 
event  that  some  parameter  needs  to  be  changed,  the  operator 
is  given  the  opportunity  to  rerm  this  portion. 

The  operator  is  instructed  to  specify  a  relative 
amplitude  (from  0  to  1.2)  that  is  less  than  the  maximum  of 
the  admittance  peak  of  interest.  The  computer  princs  rhe 
relative  amplitudes  greater  than  this  value  and  the 
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associated  frequences.  This  allows  the  operator  to  select 
the  resonance  of  interest  and  provide  the  center  frequence 
and  halfwidth  to  be  considered  in  finding  the  resonant 
frequency  and  quality  factor  .  A  comparison  is  done  to 
find  the  absolute  maximum  value  for  amplitude  and  the  half 
power  points  for  the  determination  of  a.  The  resultant 
factors  are  printed  and  the  option  to  rerun  is  again 
extended. 

This  same  subroutine  is  again  used  to  obtain  similar 
data  for  the  real  pact  of  the  iaoedance.  The  subroutine  is 
an  adaptation  of  a  part  of  a  procedure  designed  by 
Conte  [ Ref .  13 ]. 

Although  two  repeats  of  the  subroutine  take  about  10 
minutes,  the  operator  is  given  sufficient  information  on  the 
impedance  and  the  admittance  responses  of  the  transducer  to 
decide  which  type  of  measurements  will  provide  the  most 
accurate  information. 

2.  fiit a  Collection 

The  computer  has  the  frequency  synthesizer  sand  a 
frequency  far  below  resonance.  The  meters  cf  the  Dranetz 
are  zeroed  by  the  operator  and  set  for  full  scale 


deflection 


An  average  of 


ten  measurements  for  each 


Figure  3. 2  Flow  Chart  of  the  Program 
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Figure  3.2  Continued) 
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terminal  is  placed  in  memory  to  be  used  as  normalization 
factors  for  all  later  DC  output  iroltage  measurements. 

While  still  at  this  frequency  far  below  resonance, 
an  average  of  ten  aeasurements  will  be  mada,  normalized  and 
scaled  (using  operator  input  for  the  scale  factors  from  the 
Dranetz  instrument)  and  stored  as  blocked  values  of 
conductance,  susceptaoce,  resistance  and  reactance.  These 
will  be  used  later  ia  computations  after  being  assigned  to 
different  variables  depending  on  the  medium.  This  one 
measurement  is  all  that  is  necessary  for  the  electrically 
coupled  transducers  with  small  dielectric  losses.  For  the 
magnetically  coupled  transducers,  the  values  are  taken  far 
below  resonance  and  stored.  The  program  loops  back  through 
this  section  after  sending  a  message  to  the  frequency 
synthesizer  to  go  an  equal  frequency  distance  above 
resonance.  &n  average  of  the  two  values  is  stored  for  later 
use.  This  averaging  is  necessary  since  the  eddy  current  ana 
hysteresis  losses  cause  a  steady  increase  in  the  blocked 
impedance  values. 

The  decision  of  tie  operator  on  the  type  cf  data  to 
be  collected  is  entered.  The  scale  factor  from  the  Dranetz 
instrument  is  entered  and  the  oomputer  goes  to  a  subroutine 
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for  data  collection.  The  computer  draws  the  appropriate 
center  fraquency  froi  the  data  collected  during  the  search 
and  bases  the  bandwidth  for  the  sweep  on  the  quality  factor 
■easured  previously.  Fifty  frequencies  and  the  associated 
real  and  imaginary  data  are  measured  and  stored  froa  ten 
tiaes  the  half  power  bandwidth  below  center  frequency  to  two 
tiaes  the  bandwidth  below  resonance.  Two  hundred 
aeasureaents  are  talcen  in  the  vicinity  of  resonance  and 
another  fifty  in  a  region  froa  two  to  ten  times  the 
bandwidth  above  resonance.  (Note:  This  program  reduces  the 
range  of  interest  to  plus  and  minus  five  tiaes  the  bandwidth 
around  resonance  when  tests  are  done  under  load  (water) 
conditions.  The  second  program  takes  data  at  three  hundred 
equally  spaced  frequencies  within  the  bandwidth  specified  by 
the  operator.)  Dace  the  data  nave  been  collected,  another 
subroutine  is  entered  to  find  tne  minimum  and  maximum  values 
for  the  two  sets  of  data  and  their  associated  frequencies. 
This  is  accomplished  through  a  comparison  to  a  minimum  (or 
maximum)  value  until  a  smaller  (larger)  one  is  found.  A 
reassignment  is  done  and  the  comparison  continues  throughout 
the  data.  The  maximum  value  and  the  two  adjacent  values  are 
then  fitted  to  a  quadratic  equation  which  is  used  to  produce 
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a  better  estimate  of  the  maximum  and  its  associated 
frequency.  In  the  a  vent  of  an  error  in  data  collection, 
there  is  an  opportunity  to  retake  the  data  before  proceeding 
in  the  prog  ran. 

A  subroutine  provides  for  a  listing  of  the  collected 
data  (100  points)  in  the  vicinity  of  resonance.  All  of  the 
data  could  be  printed;  however,  the  data  near  resonance 
provides  the  important  information  and  printing  is  an  option 
open  to  the  operator. 

Plots  of  the  data  nay  be  made.  The  operator  may 
chose  to  make  plots  of  real  or  imaginary  amplitudes  versus 
frequency  or  the  plots  of  real  versus  imaginary  data  over 
frequency,  or  all,  or  none.  Each  plot  is  drawn  using  a 
subroutine  and  labeled  with  pertinent  scaling  and  captions. 
At  the  end  of  this  chapter.  Figures  3.3  through  3. 1U  are 
examples  of  plots  for  a  piezoelectric  transducer.  (This 
example  is  for  a  Type  100-5353-010,  Serial  No.  1213, 
designed  for  use  ia  an  active  sonobuoy.)  Figures  3.15 
through  3.22  are  examples  of  plots  for  a  magnetostricrive 
transducer.  (This  example  is  a  lype  CMT-10255,  Model  QGB, 
Serial  No.  318.) 
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3.  aaaiEaUkios. 

In  the  calculation  of  lotional  data,  the  type  of 
coupling  of  the  transducer  becoies  important, 
a.  Magnetic  Coupling 

For  a  aagnetically  coupled  transducer,  the 
notional  impedance  data  is  calculated  by  subtracting  the 
blocked  resistance  from  each  resistance  measurement. 
Blocked  inductance  at  resonance  is  calculated  using  the 
averaged  blocked  reactance  value.  The  angular  frequency 
multiplied  by  the  blocked  inductance  is  then  subtracted  from 
each  value  of  reactance.  In  the  event  admittance  data  has 
been  collected,  the  motional  conductance  is  obtained  by 
subtracting  the  blocked  conductance,  while  motional 
susceptance  may  be  obtained  by  subtracting  the  blocked 
capacitance  multiplied  by  the  angular  frequency  from  each 
value.  A  list  of  motional  data  may  be  obtained.  A  plot  may 
be  elected  and  drawn  and  labeled  with  axis  scaling  done 
using  minimum  and  maximum  values  of  the  calculated  data. 

Computations  are  done  using  the  collected  data. 
The  dynamic  electromechanical  coupling  coefficient  is 
calculated  from  the  resonant  frequencies.  A  new  quality 
factor  is  calculated  based  on  the  frequencies  of  maximum  and 
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minimum  amplitudes  of  the  imaginary  part  of  the  data 
collected.  Values  of  some  blocked  measurements,  quality 
factors,  the  coupling  coefficient,  and  resonances  are  output 
to  the  printer. 

b.  Electric  Coupling 

Por  an  electrically  coupled  transducer, 
dielectric  losses  are  small  and  it  is  acceptible  to  use  the 
values  measured  far  below  resonance  as  the  blocked 
measurements.  Notional  resistance  is  obtained  by 
subtracting  the  blocked  resistance  from  each  measurement. 
The  addition  of  the  angular  frequency  multiplied  by  the 
calculated  blocked  inductance  to  the  reactance  measurements 
will  provide  ^he  motional  reactance.  The  data  may  be  listed 
and  a  plot  made.  An  example  of  the  motional  impedance  plot 
can  be  seen  in  Figure  3.9. 

Motional  admittance  data  may  be  obtained  by 
subtracting  the  reciprocal  of  the  blocked  resistance  from 
each  conductance  value.  The  angular  frequency  multiplied  by 
the  blocked  capacitance  is  subtracted  from  each  susceptance 
measurement.  Figure  3.10  shows  a  motional  admittance  plot 
for  a  piezoelectric  transducer  in  air.  Calculations  similar 
to  those  for  magnetic  coupling  are  done  and  desired  values 


After  measurements  ire  taken  vith  the  transducer 
under  load  (in  water*  calculations  are  done  to  find  the 
efficiency  and  optimal  operating  frequency  (for  magnetic 
coupling).  This  information  is  output  to  the  printer.  An 
opportunity  to  repeat  data  collection  in  either  medium  is 
afforded  before  data  on  air  and  water  temperatures  and 
transducer  identification  numbers  are  requested  prior  to  the 
end  of  the  program.  A  sample  output  (for  the 
magnetostr ictive  transducer  Type  CUT- 10255,  Serial  No.  318, 
Model  QGB)  is  included  as  Figure  3.23. 
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Figure  3.8  El-ctrical  Input  Impedance  Plot  for  a  Piezoelectric  Transducer 


Figure  3.9  Motional  Iepedance  Plot  for  a  Piezoelectric  Ti 


Flqura  3. 11  Conductance  Spectrum  for  o  Piezoelectric  Troneducer 
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Fiqoro  3. 13  Input  Electrical  Admittance  Plot  for  a  Piezoelectric  Transducer 


R  VS.  FREQUENCY  IN  AIR 

5  Reel stance  Spectrum  for  a  Maqnetostrictive  Ti 


Fiqura  3. 16  Raoctonca 
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Electrical  Impedance  Plot  for  a  Maqnetostrictive  Tn 


VALUES  MEASURED  IN  AIR  FROM  ADMITTANCE  DATA 


BLOCKED  INDUCTANCE  *  8.86F.-4  HENRIES  AT  23956  HZ 
ELECTRICAL  QUALITY  FACTOR  =  189.73 
MECHANICAL  QUALITY  FACTOR  *  29.426 
ELECTRICAL  RESONANCE  =  23956.0379323  HZ 
MECHANICAL  RESONANCE  =  23920  HZ 
DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIENT  - 
BLOCKED  RESISTANCE  =»  33.724  OHMS 
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VALUES  FOR  MEASUREMENTS  IN  WATER  FROM  AEMITTANCE  DATA 

BLOCKED  INDUCTANCE  =•  9.01E-4  HENRIES  AT  23932.48  HZ 
ELECTRICAL  QUALITY  FACTOR  -  18.513 
MECHANICAL  QUALITY  FACTOR  »  19.688 
ELECTRICAL  RESONANCE  -  23932.4882297  HZ 
MECHANICAL  RESONANCE  -  23892  HZ 

DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIENT  *  .00338 
POTENTIAL  EFFICIENCY  -  .20398 
FREQUENCY  OF  OPTIMUM  EFFICIENCY  •  24246.61  HZ 
MECHANICAL  POWER  UTILIZATION  FACTOR  -  .14645 


Figure  3.23  Sample  Computer  Output 
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IV.  TE£T  £ND  EVALUATION 


Although  many  of  the  calculations  made  by  the  computer 
are  basic,  it  was  desireable  to  have  a  set  of  data  collected 
manually  to  ensure  accuracy  of  the  semiautoaated  system. 

A.  TESTING 

1 .  Met  hpd 

Data  were  collected  for  four  different  transducers, 
all  of  which  were  later  used  to  test  the  computerized 
system.  A  similar  equipment  set-up  was  used  with  the 
computer  replaced  by  an  operator  and  an  x-Y  plotter,  and  a 
Wavetek  pulse/function  generator  in  place  of  the 
synthesizer/function  generator. 

The  frequency  spectrum  was  swept  manually  with  plots 
made  of  the  real  and  imaginary  components  versus  frequency. 
From  these  plots,  amplitudes  and  frequency  values  were 
obtained  and  calculations  similar  to  those  done  by  the 
computer  performed. 

2.  Com  asjsgfis 

For  the  magnetostrictive  transducer  used  here,  the 
data  output  of  the  computer  corresponds  very  closely  to  that 
collected  manually.  The  blocked  inductance  was  found  to  be 


the  same  calculated  by  either  method  to  five  digits. 
Besonar.ce  frequency  variations  up  to  ten  hertz  may  be 
attributed  to  temperature  variations;  and  these  frequency 
variations  will  affect  the  coupling  coefficient.  The 
efficiencies  compared  to  within  five  percent. 

Comparisons  of  the  piezoelectric  data  ware  as 
favorable  as  that  for  the  magnet ostrictive  transducer.  The 
blocked  capacitances  varied  by  up  to  six  percent  and  there 
were  small  resonance  frequency  differences  in  some  cases. 

B.  EVALUATION 

The  resonance  frequency  variations  noted  may  be 
attributed  irainly  to  temperature  changes  from  one  experiment 
to  another.  Larger  deviations  noted  in  a  few  cases  are  due 
to  inaccuracies  in  reading  graphs  or  in  determining  a  true 
maximum  for  a  very  low  2  device  as  is  the  case  for  many 
transducers  when  placed  under  load.  For  the  examples  used 
here,  the  computer  data  are  considered  more  accurate  since 
it  measures  and  compares  data  to  five  digits,  far  better 
than  the  average  eye  can  do  with  unsealed  data.  In  the 
computerized  system,  the  data  collected  in  the  vicinity  of 
resonance  are  compared  to  find  the  absolute  maximum  value 
and  the  associated  frequency.  This  point  and  two  adjacent 
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points  are  then  used  for  a  parabolic  fit  to  find  the  precise 
frequency  for  maximum  coapooent  amplitude.  This  provides 
accuracy  in  the  resonance  frequency  far  in  excess  of  what 
aay  be  achieved  aanually  by  an  operator. 

A  'delay  time'  to  allow  stablization  of  the  systea 
(dominated  by  the  response  tiae  of  the  Dranetz  meter)  was 
determined  to  be  one  hundred  milliseconds.  In  manually 
sweeping  the  frequency  spectrum,  stablization  time  is 
difficult  to  achieve.  This  aay  result  in  the  appearance  of 
spikes  in  the  plots  which  increase  the  difficulty  of  manual 
evaluation  of  the  collected  data.  Once  this  stablization 
time  was  determined,  all  wait  times  in  the  computer  program 
were  ensured  to  be  greater  than  this  value  with  the  longest 
waiting  period  during  detemina  tion  of  the  shunt  and 
normalization  values.  The  second  longest  wait  period  occurs 
during  the  data  collection.  Here,  a  period  of  two  hundred 
fifty  milliseconds  has  been  used  in  the  vicinity  of 
resonance  to  minimize  the  chance  of  error.  The  added  minute 
for  data  collection  that  this  additional  wait  time  requires 
is  of  minor  inconvenience  compared  to  the  improvement  in 
accuracy. 
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Transducers  may  exhibit  a  moderately  high  Q  in  air 
(above  25)  .  For  these  (and  most  resonances  with  a  Q  above 
5)  the  program  with  set  frequency  ranges  for  data  taking 
works  well.  However,  these  sene  transducers  nay  have  very 
low  quality  factors  under  load  (in  water).  To  avoid  data 
review  over  a  frequency  range  beyond  the  upper  limits  of  the 
instrunents  and  perhaps  overlapping  other  resonances  of  the 
transducer,  it  is  recommended  that  the  program  allowing  for 
operator  input  of  the  frequency  range  be  utilized. 

C.  COHCLOSIONS 

It  was  noted  that  most  transducers  have  more  than  one 
resonance.  To  be  a  useful  analytical  tool,  the  operator  of 
this  system  needs  to  know  the  primary  design  frequency  of 
operation  for  the  test  transducer.  Should  the  resonance  at 
this  frequency  cause  the  admittance  or  impedance  to  exceed 
the  limits  cf  the  Dranetz  meter,  the  operator  may  opt  for 
data  taking  of  the  type  to  fall  within  limits. 

Manual  collection  of  impedance  and  admittance  data  for  a 
single  transducer  in  two  media  (air  and  water)  and 
calculations  similar  to  those  done  by  the  computerized 
system  utilizing  data  at  three  hundred  frequencies  requires 
on  the  order  of  fifteen  hours.  This  system  carries  out 


these  operations  ia  seventy- four  sinutes.  For  most 
transducers,  only  one  type  of  data  (adeittaace  or  impedance) 
is  normally  collected.  A  typical  measurement  which  produces 
results  and  the  plots  similar  to  Figures  3.15  through  3.23 
would  take  fifty-two  minutes.  (Since  each  collection  of 
data  allows  for  four  plots  and  two  lists  of  data,  the  time 
tc  run  the  program  collecting  impedance  and  admittance  data 
in  both  media  could  be  reduced  by  forty  minutes  and 
collection  of  one  type  of  data  in  both  media  by  twenty 
minutes  if  no  plots  or  lists  were  made  and  only  the  listing 
of  transducer  parameters  similar  to  Figure  3.23  desired.) 

The  major  time  saving  factor  for  this  system  would  be 
found  through  the  replacement  of  the  Dranetz  Impedance  Mater 
with  an  instrument  that  does  aot  require  as  much  operator 
interface.  However,  semi-automation  as  provided  in  this 
system  constitutes  a  major  time-saving  over  manual  means 
even  when  the  equipment  set-up  is  not  altered  in  a  major 
way. 

The  plots  produced  using  an  X-Y  plotter  and  sweep 
oscillator  enable  measurements  of  frequency  with  an  accuracy 
of  about  one  percent.  Using  the  computer,  the  precision  is 
better  than  one  ia  ten  taousand.  Similarly,  for 
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calculations  using  real  and  imaginary  components  of  the 
impedance  or  admittance,  computer  precision  is  better  than 
one  in  ten  thousand  as  opposed  to  about  one  m  ten  for 
manual  manipulations.  However,  the  measurement  accuracy  is 
at  best  plus  or  minus  two  percent  for  the  Dranetz  meter. 
The  advantage  gained  over  manual  methods  by  using  the 
computerized  system  is  due  to  the  systematic  pauses  for 
Dranetz  meter  stablization  prior  to  all  measurements. 

This  system  has  been  designed  to  analyze  any  type  of 
transducer.  Specific  adaptation  to  handle  a  particular 
transducer  type  with  known  resonance  range  can  drastically 
reduce  much  of  the  operator  interface  necessary  for 
evaluation  of  the  properties  of  the  transducer  and  reduce 
the  time  requirement  accordingly.  Adaptation  for  a 
particular  type  of  transducer  will  allow  variables  currently 
in  the  computer  program  to  be  made  constant.  Careful 
selection  cf  these  values  will  maximize  accuracy  and 
precision  in  the  desired  measurements  and  calculations. 

D.  RECOMMENDATIONS  FDR  FURTHER  IMPROVEMENTS 

In  all  of  the  preceding  discussions,  it  must  be  apparent 
that  the  "weak  link”  in  this  system  is  the  Dranetz  impedance 
meter  which  is  over  fifteen  years  old.  The  difficulties 
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with  this  device  ice  primarily  its  lack  of  computer 
controllability,  its  plus  or  miaus  two  percent  accuracy,  and 
its  high  noise  at  low  signal  labels. 

Although  it  was  unfortunately  not  available  for  this 
study,  the  Hewlett-Packard  4192A  Low  Frequency  Impedance 
Analyzer  seems  to  be  the  proper  replacement  for  the  Dranetz. 
It  is  capable  of  fill  program  control  via  the  HP-IB  and 
would  also  eliminate  the  neel  for  a  separate  frequency 
synthesizer  and  data  acquisition  system.  It  is  highly 
recommended  that  this  device  be  substituted  in  any  future 
applications  of  this  technique. 
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Two  programs  were  written  to  iapleaent  the  measurement 
procedures  and  the  calculations  to  evaluate  any  type  of 
transducer.  The  first  program  titled  MHTDR&2"  automatically 
looks  at  a  frequency  range  of  ten  times  the  bandwidth  on 
either  side  of  the  resonance  frequency.  (This  is  modified 
for  water  to  be  five  tines  the  bandwidth  on  either  side.) 
This  program  also  m toaatically  selects  the  frequency  at 
which  to  take  shunt  measurements  based  on  the  Q  factors. 
The  second  program,  titled  ’'PrCjCSH'',  allows  the  operator  to 
specify  the  bandwidth  for  analysis  and  the  frequency  at 
which  to  measure  shunt  values.  The  following  reference 
guide  is  applicable  Cor  both  programs.  Major  differences  in 
the  programs  occur  between  lines  2501  and  2650  and  between 
lines  8601  and  8635. 
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REFERENCE  GUIDE 


Line  Nos. 

1  -  60 
60  -  90 

100-149 

150 

153-176 

177-186 

198-450 

451-5  00 

525-540 

541-642 

645-657 

658-668 

670-757 

800-940 

1000-1031 

1040-1146 

1  150-1225 

1350-1410 

2501-2650 

7000-7545 

8100-81 10 
8200-8250 
8300-8380 

8400-8450 
8549- 8573 


Ds  scription 

Da clarations,  inputs,  "bookkeeping” 

Admittance  Spectrum  (fcsq.  6  Q)  using 
subroutine  7000 

Iipedange  Soectrum  (freq.  &  Q)  using 
subroutine  7300 

*jptaias  blacked  values  using  subroutine 

Impedance  data  collection  (option) 

Admittance  data  collection  (option) 

Listing  and  plots  of  data  using 
subroutines 

Notional  data  calc,  for  magneto,  (for 
’!•  data) 

Notional  data  calc,  for  elec,  (for 
’?•  data) 

List  of  motional  data  and  plot 

Notional  data  calc,  for  elec,  (for  ' Z' 
data) 

Notional  data  calc,  for  magneto,  (for 
•Z  •  data) 

List  of  motional  data  and  plot 

Calculations  for  elec,  coupling  in  air 

Calculations  for  magn.  coupling  in  air 

Calculations  for  elec,  coupling  in  water 

Calculations  for  magn.  coupling  in  water 

Options  for  rerun  or  change  of  medium 

Subroutine  to  determine  blocked  values 

Subroutine  for  search  and  spectrum 
analysis 

Subroutine  to  print  data  list 

Subroutine  to  plot  real  data  vs.  freq. 

Subroutine  to  olot  imaginary  data  vs. 
creq. 

Subroutine  to  plot  ’circles' 

Subroutine  to  obtain  normalization 
factors 
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I 


8600-8699 

8700-8799 

9000-9350 


Subroutine  to  collect  data 
Subroutine  to  determine  max/iin  values 
Subroutines  for  printouts 


1  REM  **HYDRA2** 

2  OPTION  BASE  I 

3  SHORT  GO  ,vl,Q8  ,0i9  ,w  ,£  ,E2  ,E3,E5  ,R8  ,R7,X2,A2,X5,X6,X?,X8 
,G7,G8  ,£2  ,£  3 ,£5  ,£6  ,E7  ,£8  , K  ,  Kl  ,P  5  ,P6  ,w 

7  SHORT  A  ,F2  ,F3  ,F4  ,F5  ,H3  ,H4  ,J1  ,02  ,U  3,  J4  ,R  ,X  ,G  ,E 
6  INTEGER  M  ,  Y  ,  Y 1  ,1  ,Y2  ,J  ,22  ,fi  1  ,H2  #L  2  ,S  I 

10  DIM  fc$[2]  ,A§  [  2  ]  ,T$  ( I J  ,D(.3CG,3)  ,A1$[4J  ,D$(30Cj  ,L$[20J  , 
F.Sll]  ,15(1]  ,S5(1U]  ,C$  [  5]  ,  G  $.  [  2  3  ]  ,H${23] 

11  REAL  S  ,N  ,N8  ,N9  ,R0  ,R9  ,X0  ,XS  ,L  ,C0  ,L0  ,0  ,P1  ,H6  ,HS  ,F,W0  ,GG 
,G9  ,BG  ,E9  ,N(  10)  ,M(IG)  ,F0,FL,F6  ,F6  ,F9 

14  REAL  S6  ,  E4 ,E6  ,L5  ,D  ,D1  ,D2  ,D  3  ,D4 

15  Xb  ,D  1  ,D2  ,L  3  ,C4=0 

16  t,E2  ,E3,E5,R6,R7,X2,X3,X5=0 

17  2.3  ,A,F2  ,F3  ,  F*=  ,F5  ,B6  ,  B7  ,  £6=0 
IS  SS  ,£4  ,ES  ,M  ,  Y  ,Yl  ,1  ,Y  2  ,J  ,  22*0 

19  K  ,Kl  ,F5  ,F6  ,W  ,00  ,01  ,0'S  #09  ,E=0 

20  CO  ,LC,C  ,P1  ,E,K8  ,H9  ,P  ,W0  ,S*0 

21  N  ,NS  ,N9  ,R0  ,R9  ,X0  ,XS  ,G0  ,G9=0 

22  £C,89,Hl,H2  ,F0  ,Fl  ,F6  #FS  ,F9=G 

25  DISP  "THIS  PROGRAM  AUTOMATICALLY  SELECTS  HALFWIDTH  FC 
E  DATA  COLLECTION.  " 

26  LISP  "THIS  LOOKS  AT  10  X  EW  IN  AIR  &  5  X  EW  IN  WATER. 

II 

27  LISP  "TO  INPUT  YOUR  OWN  BANDWIDTH  USE  'PICKBW'  PECGKA 
M  INSTEAD.  (HIT  "CCNT')" 

29  PAUSE 

30  CLEAR  (t  BEEP 

.35  DISP  “ENTER  TYPE  OF  TRANSDUCER  TO  BE  MEASURED ?  (  V '  FC 
h.  MAGNETIC  CCOFLING  OR  ' 'E  '  FCR  ELEC.)" 

42  INPUT  TS3  DISP  "TC  GET  A  COMPLETE  SET  OF  DATA  YOU  NEE 
C  MEASUREMENTS  IN  EC^S  AIR  AND  WATER .  " 

45  DISP  "DC  AIR  FIR!  aK  WHAT  MEDIUM  ARE  YOU  OPERATING?  ( 
ENTER  'l  '  FOR  AIR  OF  '2  '  FCR  WATER)" 

<-.6  INPUT  M(s  CLEAR  §  BEEP 
47  IF  »2  THEN  49 
S  C  S  =  "  A I R  "  1?  Z2=  10  (£  GOTO  51 
49  C$=" WATER"  U  Z2  =  5 
51  PRINTER  IS  2  is  PLOTTER  IS  1 

55  LISP  "SET  THE  CRANETZ  FCR  ADMITTANCE  .SET  ON  LCwEST  SC 
ALE  THAT  wCN  'T  PEArv  CURING  RUN.iCONT)" 

6  C  PAUSE 

o 5  PRINT  iS  PRINT  s  PRINT 

70  PRINT  "ADMITTANCE  IN  “;CS 

75  GOSUE  7000 

60  IF  M=2  THEN  SO 

65  FU*F6  @  wG=0  <s)  GOTO  ICO 

90  F1=F6  is  ul-w 

100  GCLEAR  y  CLEAR  BEEP 

105  DISP  "SET  DFANETZ  FCR  IMPEDANCE- ( SCALE  NEEDED) -HIT  C 
ONT  " 


110  PAUSE 

115  PRINT  £  PRINT  §  PRINT 
120  PRINT  "INPEDANCE  IN  " ;C$ 

125  CCSUB  7000 

130  IF  fr' =2  THEN  140 

135  F8*F6  @  <s  GOTO  141 

140  F9*F6  (? 

141  CLEAR  <?  EEEP 

142  PRINTER  IS  701,76  &  PRINT  USING  144 
14  4  IMAGE  2/ 

145  PRINT  “F(Y-AIR)*  "  ;  I  NT  (F  C)  ,  F  (  Y.-WA1 )  =  "  ;  INI  ( F  3 )  ,  ,"F 
(Z.-AIR)  =  ”  ;  IN.1  ( F8 )  ,  F  ( 2 -WAT)  =  " ;  INI  (F  5 ) 

146  PRINT  " w ( Y-AI K ) *  " ; GO ,  w ( V- WATER ) =  " J vl  ,  ,  "W ( Z-AIR) • 
";CS  ,  ,"C  (Z-kvATER)  =  ";wS 

147  PRINTER  IS  2 

148  CLEAR  y  DISP  "THESE  ARE  THE  INITIAL  ESTIMATES . (HIT  " 
CCN'l"  rtHEN  READY  TO  PROCEED) “ 

149  PAUSE 

150  GCSUE  2501 

152  CLEAR  <y  EEEP 

153  DISP  "TO  COLLECT  ADMITTANCE  DATA  (FOR  ELEC.  CCOFLING 
)  ENTER  'l  ' .  ENTER  '2'  FOR  IMP  ELAKCE" 

154  INPUT  L2<s  PRINTER  IS  2 

155  IF  L2=2  THEN  158 

156  R$="G"  d  (i  S  $=" ADMITTANCE"  H$="CCNEUCTANCE 

(MICROMHOS)"  £  G$*"SUSCEPTANCE  (MICRCMHCS)" 

157  GOTO  159 

158  R$  =  "R "  t  I$«"X"  fc  S$="  IMPEDANCE"  1$  G$=" RCACTANC E  (CH 
MS)"  6  H$»"RESISTANCE  (OHMS)" 

159  IF  L2-1  THEN  176 

160  CLEAR  y  DISP  "SET  ON  Z  AND  ENTER  SCALE  FACTCF  TO  COL 
LECT  DATA." 

161  INPUT  S6@  PRINT  SS;"  IN  " ;  CS  C°  PRINT  "SCALE  =";So;" 
OHMS  " 

162  CLEAR  y  DISP  "I  AM  FORKING" 

164  GCSUE  8600 

165  IF  Y=2  THEN  152 

166  IF  M*2  THEN  165 

167  R&=£  l?  X3=E2  <s  X  o  =E  3  y  F  4=E4 

168  X6=E  5  y  F5  =  E6  <s  F6  =  E1  y  GCTC  198 
lo9  P7  =  E  is  X2  =  £2  <s  X7  =  E3  y  F4=E4 

170  X5=£5  y  F5=E6  b  FS=El  y  GOTO  19b 

176  CLEAR  y  SEEP 

177  DISP  "SET  CN  Y  &  ENTER  SCALE  FACICR  IN  MICRCMHCS  TC 
COLLECT  DATA." 

178  INPUT  Sby  PRINT  SS;"  IN  ";C$  y  PRINT  "SCALE  = " ; S 8 ; " 
MlCRCMhOS" 

i  79  GCSUE  8600 
18U  IF  Y  =  2  THEN  152 
lbl  IF  M=2  THEN  185 
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162  G6«£  @  B3«E  2  @  £8*E3  @  F2*E4 
363  36*E 5  @  F3»E6  te  F0*E1  (i  GCTC  198 
165  G7«E  €  B2-E2  g  B7=E3 
186  F2»£4  @  B5*E5  @  F3*E6  §  Fl»El 
196  PRINT  USING  199 

199  IMAGE  2/ 

200  OISP  "DO  YOU  DESIRE  A  LIST  OF  COLLECTED  DATA  NEAR  RE 
SONANCE?  (100  POINTS )  ( 1*YES  ,2=NC ) " 

201  INPUT  Y 

202  IF  Y*2  THEN  300 
204  IF  L2=2  THEN  240 

210  PRINT  "ADMITTANCE  DATA  FOR  “;C$;"  IN  MICROMHOS*'  &  PR 
INT  USING  215 

215  IMAGE  2/ FREQUENCY"  ,  4X  ,  "REAL"  ,6X  , "IMAGINARY" 

216  GOTO  265 

240  PRINT  "IMPEDANCE  DATA  FOR  “  ;C$;  "  IN  OHMS'*  (<?  PRINT  US 
ING  215 

265  GOSUB  8100 

300  CLEAR  &  DISP  "DC  YOU  DESIRE  A  PLOT  OF  DATA? ( 1=Y£S ,2= 
NO)  " 

30  5  INPUT  Y@  GCLEAR  CLEAR 
315  IF  Y=2  THEN  450 

.320  DISP  "ENTER:  1=  G/R  VS .  FREQ;  2=  B/X  VS.  FREQ;  3=  B 
VS.  G/X  VS.  R;  4*  END  PLOTTING  LOOP." 

325  INPUT  Y 
3 3U  IF  Y-2  THEN  365 

325  IF  Y=3  THEN  425 

340  IF  Y=4  THEN  450 

345  GOSUB  6200 

356  LDIF  0  t  PEN  1  fe  PENUP 

357  MOVE  C(l,l)+. 2*  (O  (  30  0 , 1)  -D  (1,1)  )  ,  -  ( E*  .  6 ) 

36 C  LABEL  R$;M  VS.  FREQUENCY  IN  ";C$ 

361  GOTO  320 
365  GOSUB  8300 

397  LOIR  C  «  PEN  1  <§  PENUP 

398  MOVE  D  ( 1 , 1)  +  .  2*  (D(.  300,1) -D(  1,1)  )  ,-L-.5*L 
400  LABEL  15;"  VS.  FREQUENCY  IN  ” ;C5 

4 C 1  GOTO  320 

42  5  GCLEAR  fe  CLEAR  t?  CCSUB  6400 
431  PEN  1 

442  LDIR  C  <§  PEN0P 

443  HOVE  -  ( . 5 *L)  ,-L-.2*L 

445  LABEL  "INPUT  ELECTRICAL  ";S$;"  PLOT  ( "  ;  I  5 ; " vS . " ; R$; " 
)  IN  ";C$ 

446  GOTO  320 

450  IF  L2*2  THEN  645 

451  IF  T$="£"  THEN  525 

452  GCLEAR  b>  CLEAR  @  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 
A" 

465  C0*B*.G00C01/(2*PI*F6) 
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466  L0»1/((2*PI*F6)"2*CU) 

. 470  FOR  I»1  TO  300 
475  D(I,2)«D(I  ,2).-G 
480  D(I  ,3)«D(I  ,2).-D(I  ,1)*B/F6 
485  NEXT  I 
490  GO SUB  8700 
500  GOTO  541 

525  GCLEAR  @  CLEAR  @  DISP  "I  AW  CALCULATING  MOTIONAL  DAT 

A" 

532  C0*B* . 00000 1/ ( 2*PI *F ) 

535  FOR  1*1  TO  300 

536  D  ( I  ,  2 ) *D ( I ,2 ).-G 

537  D(I,3)*D<I  ,3).-D<I,l)*B/F 

538  NEXT  I 

540  GOSUB  6700 

541  CLEAR  &  GCLEAR  @  BEEP 

542  DISP  "DC  YOU  DESIRE  A  LIST  OF  MOTIONAL  DATA  NEAR  PES 
ONANCE? ( 1=YES ,  2*NC)“ 

543  INPUT  Y 

545  IF  Y*2  THEN  604 

575  PRINT  "MOTIONAL  ADMITTANCE  DATA  IN  "  ;C$  @  PRINT  USIN 
G  580 

580  IMAGE  2/, " FREQUENCY" ,  4X  , "REAL"  ,6X  ,  "  IRAGINAPY  " 

595  GCLEAR  €  CLEAR 

600  GOSUB  8100 

6C4  CLEAR  €  BEEP  @  BEEP 

605  DISP  "WANT  A  PLOT? (  1*YES  , 2*NC ) " 

610  INPUT  Y2 

b!5  IF  Y2=2  THEN  755 

t>  16  PLOTTER  IS  705  @  PEN  1 

620  GOSUB  8400 

628  LDIR  0  g  PEN  UP  @  PEN  1 

629  MOVE  -  (. 5*L)  ,-L-.2*L 

630  LABEL  "MOTIONAL  ";S$;"  PLOT  FOR  ";CS 
642  GOTO  755 

645  CLEAR  fc  GCLEAR  {?  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 
A" 

646  IF  TS*"K"  THEN  656 

647  FOR  1*1  TO  300 

648  D( I  ,2) *D( I ,2) -R 

649  CG»1/(2*PI*F*X) 

650  L0* 1/ ( (2*PI*F6)~2*CQ) 

651  D(I  ,3)*D(I  ,3).-2*PI*D(I  ,1)*L0 
o  52  NEXT  I 

653  GOSUB  8700 

654  CO-AES(CU)  <?  LG*ABS(L0)  «  BEEP 

655  DISP  "DO  YOU  WANT  A  LIST  OF  THE  DATA  NEAR  RESONANCE? 
( 1*YES  ,2*NO) " 

656  INPUT  Y 

657  GOTO  670 


658  L0-X/(2*PI*F6) 

662  FOR  1*1  TO  .200 

663  C(I  ,2)*D(I  ,2) -R 

664  D ( I , 3 ) *D ( I  ,  3).-C  ( I  ,  1)  *2*PI*LG 

665  NEXT  I 

666  GCSUB  8700 

667  BEEP  @  DISP  "  DC  YOU  WANT  A  LIST  CF  DATA  NEAR  RESCNA 
NCE?(1*YES,2*NQ)“ 

668  INPUT  Y 

670  IF  Y*2  THEN  706 

665  PRINT  "MOTIONAL  IMPEDANCE  DATA  IN  ";C$  £  PRINT  USING 
666 

666  IMAGE  2/ FREQUENCY" , 4X  ,MR£AL"  ,6X , "IMAGINARY" 

7  00  GCLEAR  £  CLEAR 

705  GCSUB  6100 

706  CLEAR  £  BEEP  £  BEEP 

710  DISP  "WANT  A  PLOT? (  1=YES  , 2=NQ ) " 

715  INPUT  Y2 

72C  IF  Y2*2  THEN  755 

721  PLOTTER  IS  705  £  PEN  1 

725  GCSUB  6400 

737  LDIR  0  £  PENUP 

736  MOVE  -  (.5*L)  ,-L-.2*L 

740  LABEL  "MOTIONAL  " ;S$;”  PLOT  IN  ";C$ 

755  GCLEAR  £  CLEAR 

756  PRINT  USING  757 

757  It  AGE  3/ 

600  DISP  "I  AM  DOING  CALCULATIONS  FOR  YOU.  PLEASE  EE  PAT 
I  ENT." 

60  3  IF  xV*2  THEN  1040 
olO  IF  T5*"tf"  THEN  100U 
825  IF  L2=2  THEN  860 
84C  w0=El/AES  (E4-E6) 
o45  CL*A6S  (E2/ (2*PI*FG)  ) 

655  GOTO  665 

66  0  y 6=E l/ABS (E 4-E6 ) 

665  C1*A£S( 1/ (2*PI*FS*E2) ) 

665  K= 1- (FC/Fb ) * 2 

655  Ki*Cl/(ABS(CC)+ClJ 

S:C  GCSUB  9GGG 

921  PRINTER  IS  2 

940  GOTO  1250 

lGOG  K* 1- ( F6/F0 )  "2 

LuuS  IP  L2*2  THEN  1g25 

1010  y G*£ l/AES (£ 4-E6 ) 

lul5  GOTO  lu30 

1025  yb*£l/A£S (E4-E6) 

1g30  GCSUL  820g 

lo 3 1  PRINTER  IS  2  GOTO  1350 

lo40  IF  THEN  1150 
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1C41  D ISP  "IF  YOU  HAVE  DATA  FOR  AIR,  ENTER  1;  ELSE  2" 

104 2  INPUT  Y 

1043  RAD 

1044  IF  L2»l  THEN  1IIC 
IC45  IF  Y«2  THEN  1062 

1060  P6*D4*  (D3-D4)/ (D3*R7) 

1061  wOFl 

1062  wl“E  1/ASS  (E  4-E6)  @  K* 1-(F  1/F  9)  *2 
1070  GOSUB  S10C 

1076  PRINTER  IS  2  £  GOTO  1350 
1110  IF  Y*2  THEN  1135 
1126  NC-F1 

1130  F6«C2*  (D1-D2)/(D1*C7) 

1135  U1*E 1/AES {£  4-E6) 

1140  K*l- (Fl/FS) *2 

1145  GCSUB  9100 

1146  PRINTER  IS  2  @  GOTO  1350 

1150  RAD 

1151  DISP  "IF  YOU  HAVE  DATA  FOR  AIR,  ENTER  1;  ELSE  2" 

1152  INPUT  Y 

1153  IF  L2* 1  THEN  1160 

1154  IF  Y*2  THEN  1162 

1155  B».5*ACS( (R7-R9)/D4) 

1156  H6*SIN<&)*2  y  HS=COS(S)~2 

1157  P5®(  ( 1+D  3/R0*H  9)  5- (  L-D  3/R0*H8  )  *.5)/(  (1+D3/R0*H9)  “ 

. 5+ ( 1-C  3/R0*H8)  A  .5) 

1159  P=D4*SIN  (2*B)/ (4*R9*QS)  @  aG®  F9*  (2*P  +  (P  *2+1)  "  .  5) 

1160  P6®L4* (D3-D4)/ (D3*R7) 

1162  uS*El/ABS(E4-E6J  y  K»l-(f  S/Fi)  *2 

1163  PRINT  "E® " , "K= " ; K , " hi® " ; R 1 , " H6=" ;H6 ."HS®" ;H5 

1164  GGSUE  9300 

1166  PRINTER  IS  2  t  GOTO  1250 
1160  RAD 

1162  IF  Y *2  THEN  1210 
1164  £®.5*ACS(  (G7-GS.)/D2) 

1186  H8*SIN  ( B )  *  2  @  h9=CCS  (5)  "2 

1 188  P  5*  (  (1+D1/GC*HS)  5-(  1-0  1/G0*H6  )  *.  5)/  (  (1+Cl/GU*H  5)  * 

. 5+ ( 1-Dl/G0*H8) *  .5) 

1190  P*D2*S1N  (2*E)/ (4*G9*^1)  is  *v0®F  1*  (2*P  +  (P  '2+1)  *, .  5) 
1195  P6«C2* (C1-C2)/ (C1*C7) 

1210  W5«E1/AES  (E4-E6)  y  K  =  l-{  F  9/F 1 )  “2 

1211  PRINT  "E=" ;£  , "K®”; R  R  1="  ;K1,"H6="  ;H6  ,MH9®"  ;HS 
1215  GCSUB  9300 

1225  PRINTER  IS  2  g  CLEAR  y  EEEP 

1350  DISP  “DC  YOU  DESIRE  ANOTHER  IN  THIS  NEDIUM? (TO  GET 
THE  OTHER  TYPE  DATA? )( 1=Y ES  ,2®NC ) " 

1355  INPUT  Z1 
1357  PLOTTER  IS  1 

1360  IF  21*1  THEN  152 

1361  PRINT  USING  1362 
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1262  IMAGE  3/ 

1365  DISP  "DC  YCU  DESIRE  A  RUN  IN  ANOTHER  MEDIUM? ( 1=YES , 
2*NO)  " 

1370  INPUT  Z1 
1375  IF  Zi*l  THEN  45 

1360  DISP  "JUST  FOR  THE  RECCED ,  INPUT  WATER  TEMP,  AIR  TE 
MP,  TRANSDUCER  SER  NO,  MODEL  NC.,M 
1381  DISP  "AND  TYPE  (E  OR  M) “ 

L265  DISP  "INPUT  'o'  IF  IKFC  IS  UNKNGWN" 

1350  INPUT  K  ,A  ,S  ,N  ,L$ 

1351  PRINTER  IS  701,76 

1400  PRINT  "WATER  TEMP*" ; n  ,  ,"AIR  TEMP*" ;A  ,  SER.  NG.  =  ";S 
,  ," MCDEL  NC.*H  ;N  ,  ."TYPE  ";L$ 

1410  GCTO  5956 
2501  REM  *ShUNT* 

250  3  CLEAR  <s  SEEP  t§  DISP  "WE  ARE  FINDING  SHUNT  VALUES  FC 
R  G/E  CR  R/X." 

2510  A$="FR"  &  BR="hZ"  &  Sl=l 

2511  IF  T$*"M"  THEN  2520 

2512  IF  M*2  THEN  2516 

2512  G*MAX(U0,GS)  @  F*F0- (Z 2-2)  *F0/G 

2514  IF  F>200  THEN  2535 

2515  F=20G+(FG-2*FG/w-2C0)/50  t?  GOTO  2535 

2516  W*MAX  (ul,wSi)  f§  F=F  ]- (Z 2-2 )  *F  1/w 

2517  IF  F > 2 0 0  THEN  2535 

251c  F*20G+(F1-2*F1/w-2GC)/50  $  GOTO  2535 

2520  IF  M*2  THEN  2524 

2521  W=MAX(UC,G8)  §  F=Fd-  (4  2-2)  *fo/v 

2522  IF  F >20 0  THEN  2535 

2523  F  *20  0+  (F8.-2*  Fb/C-2  00 ) /  50  @  GCTC  2535 

2524  v®MA*  (Gl  »u9.)  £  F=F  5-  (Z  2-2}  *FS/v 

2525  IF  F>2C0  THEN  2535 

2527  F*2CG+(F9-2*fS/c~2CC)/ 50 
2525  D$=vAL$ ( F ) 

2527  OUTPUT  717  ;A$,DS,E$ 

2536  WAIT  10C0 

2540  DISP  "SET  DRANETZ  FREv .  SCALE  TC  COVER  t'REv.  CN  SiN 
TKESIZER  .ZERC  METERS  .SET  ON  FS.(CCNT)" 

2541  PAuSE 

25<t2  DISP  "HERE  aE  GET  THE  NORMALIZATION  FACTORS." 

2542  Output  709  ;"vtl" 

25<*6  GCSUB  6549 

2550  G,o,K,X=U 

2551  A6*"FR"  fc  6 $  = " H Z "  fc  D$*vALS(F) 

2553  OUTPUT  717  ;A $  ,D$  ,B$ 

2561  CLEAR  (g  DISP  "SET  DRANETZ  ON  i  ,S£T  SCALE  , SET  CL  NCR 
M  FOP  FILTER  +  rHASE.PLlG  IN  TRANSDUCER." 

2  562  SEEP  $  DISP  "EN'iER  SCALE  FACTOR  IN  MICROPROS .  SAi-iL  S 
CALL  AS  t-CR  DATA .  («E  GET  SnUNT  G/n)" 

2564  INPUT  S6(s  CLEAR  is  DISP  "I  AM  aCRKING  TC  GET  SHUNT  v 


ALOES"  b  E££P 

2565  FOR  1=1  TO  10 

2566  OUTPUT  717  ;A$,C8,E$ 

2567  *AIT  1000 

2566  OUTPUT  709  ;"AC3" 

2569  ENTER  709  ;  M(I) 

2570  G*G+N(L) 

2571  OUTPUT  709  ;"AC4" 

2572  ENTER  70S  ;  N{I) 

2573  B=B+N(L) 

2574  NEXT  I 

2575  G*G/(I-1)  g  E=B/{I-j) 

2576  G*G*S8/Nb  b  £=B*£S/NS  fc  £££P 

2576  DISP  " £ET  CRANETZ  ON  Z;SET  SCALE  FOR  RA X  RESPONSE.!*, 
GVE  TPAN’S  .  INPUT.  ENTER  SCALE  FACTOR" 

2579  INPUT  Set  CLEAR  @  DISP  "I  AF  NO F KING  IC  GET  SFUNT  v 
ALUES" 

2580  FCR  1=1  TO  10 

2561  OUTPUT  717  ;A £,£■£,£$ 

2582  rtAIT  1CC0 

2563  OUTPUT  709  ;"AC3" 

2584  ENTER  70S  ;  K(I) 

2565  R=R+1>i  ( L) 

25to 6  OUTPUT  709  ;"AC4;: 

2567  ENTER  70S  ;  N(I) 

2566  X*N(I)+X 
256 S  NEXT  I 

25SG  F=R/(I-1)  b  X=X/ ( I -  3 ) 

2591  R=R*S6/N3  (?  X=X*S6/NS 
25S2  IF  T$*"M"  TfcEN  2olG 

2555  IF  I*.=  1  THEN  26ue 

2556  RS=R  b  XS=X  b  G9=C  b  £ S=B  y  GOTO  2601 
2 feOO  GG=C  b  E0=E  b  R0=R  y  Xo  =  X 

2601  PRINT  "Siiv.u  '/ALOES " 


26C2  PRINT  USING  2603 
260  3  I  RAGE  1/ 

2605  PRINT  "GO*  "  ;GC  ,  ,  "£C  =  ";B0  ,  f  0=  "  ;F0  ,  XG=" 
;G  S,  ,“L  S="  ;  OS  ,  R  S=  "  ;  R  S  ,  , "  >.S=  "  ;  XS 


XG  ,  , " G S * " 


2.6G6  PRINT  y  GOTO  26  30 
261C  IF  £1=2  THEN  2620 

2611  J  1*G  b  02=E  i?  o  3=R  y  J4  =  X 

2612  IF  R=2  THEN  2615 

2613  F.*F6+ (Z2-2)  *Ffa/vi  y  Sl=2  b  GCTC  2550 
2615  F*F9+ (Z2-2) *F9/W  b  Sl=2  b  GCTC  2550 

2620  G* .  5*  (G+G  1 )  y  6=.5*{6+J2)  b  R*.5*(F+J3)  y  X=.5*(u4  + 

X)  b  GCTC  25S5 

2630  GCLEAR  b  CLEAR  y  EEEP 

2635  DISP  "ENTER  'l'  IF  ALL  IS  vvELL;  '2'  IF  *CG  NEED  A  R 


EPEAT  ." 

2640  INPUT  ¥ 
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2645  IF  Y*2  THEN  2501 
2  650  RETURN 

7000  REM  **P-Q** 

7001  GCLEAR  @  CLEAR 

70C  5  DISP  “ENTER  LOWER,  UPPER  FREO •  FOR  SWEEP. IF  PESCNAN 
CE  UNKNOWN  USE  10000,200000.  (F1,F2)“ 

7010  INPUT  F4,F2 

7011  DISP  "ENTER  RMS  VOLTAGE  NEEDED.  ( 30 CO  MV  FOR  LOW  F RE 
C.  AS  LOW  AS  2000  MV  FOR  ttlGH)" 

7012  INPUT  A£  A 1$*VAL$ (A) 

7020  CLEAB  fe  BEEP 

7025  A$*"AK"  «a  B$*"MR“ 

7035  OOTPUT  717  ;A$,A1$,3$ 

7040  OOTPUT  709  {“AClvTS'' 

7045  F3*CEIL(  (F2-F4)/3GU) 

70b0  CLEAR  <?  DISP  ”  I  AM  WORKING  TC  GET  THE  SPECIF UM  DATA 
FOR  YOU , THEN  WILL  SfcOW  &  MAKE  A  PLOT" 

7U65  FOR  1*1  TO  300 
7070  D (I , 1 ) *F 4+1 *F 3 
7075  A  $  *  "  H  2  "  £  BS-"FR" 

706  3  D$*VAL$  ( C  (  I  ,1)  ) 

7065  OUTPUT  717  ;BS,0$,A$ 

7C9C  OUTPUT  709  ;"AC3VT3“ 

7095  *AIT  100 

7100  ENTER  709  ;  L(I,2) 

7105  NEXT  I 

7110  F5*F4-.  1*  (F2-F4)  (°  F 0*F 2+ .  1*  (F 2-F 4 ) 

7115  Fb*F2+> .  1*  (F2-F4) 

7115  PLOTTER  IS  1  e  GCLEAR  &  CLEAR 

7125  SCALE  F5  ,F  6  ,~ .  1  ,1 . 2 

7130  XAXIS  C,2GG0,F<i,F2 

7135  YAaIS  F4,  .1  ,-,l  ,1.2 

7140  FOR  1*1  TC  300 

7145  PEN UP 

7150  PLOT  D ( I  ,  1 )  , D ( I  ,  2 ) 

7155  NEXT  l 

7156  GRAPH  (y  COPY 

7160  DISP  "NEED  ANOTHER  RUN?  {  1*YES  ,2*NC  )  " 

7161  INPUT  i 

7162  IF  i=l  THEN  7005 

7163  CLEAR  <y  GCLEAR  y  SEEP 

7165  CISP  "ENTER  DECISION  POINT  FOR  AMPLITUDE  (0  10  1.2) 
.(LEES  THAN  THE  FAX  DISPLA YEL )  ." 

7170  INPUT  S(y  CLEAR  «>  u££P 

7160  PRINT  "AMPLITUDE  IN  VOLTS";"  FREQUENCY" 

7101  PRINT 

7150  FOR  1*1  TC  300 

7155  IF  D (I ,2) <S  THEN  7210 

72ou  PRINT  USING  7205  ;  D(I,2),D(I,1) 

72  u  5  IMAGE  IX  , D  .  LCDD DDE D  ,  1  OX  .DDDDDD  .DD 

102 


V.-  I  ■  Mu'**  «Vr-  ■  ■»  - 


m 


72 1C  NEXT  I 

7215  PRINT  USING  7220 

7220  I  RAGE  3/ 

7225  CLEAR  fe  BEEP 

7226  DISP  "IF  *00  ARE  READY  TC  CONTINUE  ENTER  '1';  OTHER 
nlSE  ENTER  '2  ' 

7227  INPUT  Y 

7228  IF  Y*2  THEN  7160 

7235  DISP  "wHAT  IS  CENTER  FREQ,  HALFWID1H  TC  CONSIDER?" 

7236  INPUT  F,Nb  CLEAR  b  BEEP 

7240  DISP  "I  AM  WORKING  TO  GET  GOOD  FREQ  AND  Q  DATA’’ 

7241  FOR  1*1  TO  300 

7242  D(I,1)=F.-N+I*N/15C 

7243  A$*”HZ H  $  6$*hFR" 

7245  D$*VAL$ ( C ll ,1) j 
7255  OUTPUT  717  ;B$,D$,A$ 

7260  OUTPUT  709  ;"AC3VT3" 

7265  WAIT  100 

7270  ENTER  70S  ;  D(I,2) 

7275  NEXT  I 

7285  CLEAR  g  GCLEAR  g  DISP  "I  AM  FINDING  THE  AoSCLOTE  M A 
X  AND  FREQ  UPPER  AND  LOWER" 

72SC  fi  1 ,34*0  ig  H  ,H  1  ,H2  ,H  3  ,H  4=  1 

72S5  E2,B3*50 

7300  FOR  1*2  TC  20U 

7310  TF  D  (I  ,2) <D  |H  ,2)  THEN  7  340 

7315  IF  D  ( I  ,2) *D  (H  ,2 )  THEN  7  325 

7325  A6*D  (1,2)  g  F6*C(I,1)  g  H*I 

7340  NEXT  I 

7345  A7=A6/2 

7350  FOR  1*1  TC  H 

7355  IF  A7*C (1,2)  THEN  7405 

7360  IF  A7<D( I  ,2)  THEN  7365 

7  365  IF  D  (I  ,2) <£ 1  THEN  7425 

7  370  B  1*D  (1,2)  0  Hl  =  I  b  GOTO  7 4, 2  5 

7365  IF  C  ( 1 ,2)  >£2  THEN  7400 

7 3SC  B2*D(  I  ,2)  <g  H2*I 

7400  GOTO  7425 

7405  B 1  ,B2*D  (1,2)  g  Hl,H2=I  e  F7*C(n2,l)  0  GCTC  7440 
7425  NEXT  I 

7430  X*(A7-E1)/(E2-B1) 

743  5  F7*X*  (D(H2  ,l)-C(Hl  ,1)  )+C  (til.l) 

7440  FCR  I =H  TC  300 
7445  IF  A7*C( I  ,2)  THEN  74S5 
7450  IF  A7>D( I  ,2)  THEN  747  5 
7455  IF  D  ( I , 2)  >B 3  THEN  7470 
7460  B3*D(I,2)  g  H3*I 
7470  GCTC  7515 

7475  IF  D(I,2)<*£4  THEN  74S0 
746o  B4»C  (1,2)  It  H4*I 
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7490  GOTO  7515 

7495  B3,B4*D(I,2)  @  H3,H4*I 

7505  L5*D(H3,1) 

7510  GOTO  7530 
7515  NEXT  I 

7520  X*(A7-B4)/(B3-B4) 

7525  L5*'(X*(D<H4,1).-D(H3,1).))+D(H4,1) 

7530  C*F6/(L5-F7) 

7535  PRINT  "CENTER  FREQ  IS  ";F6 

7536  PRINT  "Q  IS  ";Q 

7536  DISP  "IF  YOU  ARE  READY  TO  PROCEED,  ENTER'l'.  TO  RE 
RUN  FOR  BETTER  VALUES ,  ENTER  '2'." 

7539  INPUT  Y 

7540  IF  Y*2  THEN  7235 
7545  RETURN 

8100  REM  **DATA  LIST* * 

6101  PRINTER  IS  2 

6102  FOR  1*99  TC  199 

6103  PRINT  USING  8105  ;  D(I,1)  ,D  ( I  ,  2 )  ,  D  ( I  ,  3  ) 

8105  IMAGE  CDDDDD  .D,2X  .D.DDDE  ,2X,D.DDDE 

8106  NEXT  I 

8107  PRINT  USING  8106 
810b  IMAGE  3/ 

8110  RETURN 

8200  REM  **PLGT  RE** 

6201  GCLEAR  is  CLEAR 

8202  H8*D(1,1)-.15*  (D(  300,1 ) -D  (1,1)  ) 

8203  H9*D ( 300 , 1)  +  • 1* (D (  300 , 1) -C ( 1 , 1 ) ) 

62C 5  PLOTTER  IS  70  5  is  PEN  1 

8209  SCALE  H8  ,H9  , - ( . 8*E )  , E  + . 25 *E 

6210  XAXIS  0,  (D(3GG,  1 ) -C  ( 1 ,  1))/1G,  C(l,  1)  ,0(300,1) 

8211  YAXIS  D(  1,1)  ,t/3  ,-(.  6*E)  ,E  +  .2*E 
6212  PENUP 

6215  FOR  1*1  TO  300 
8216  PLOT  D(I,1) ,0(1,2) 

6217  NEXT  I 

6216  PENUP  @  DEG  @  LOIR  C,SIN(9U) 

6221  FOR  L1=C  (1,1)  TO  0(  300,1)  STEF  (D  (  300 , 1 ) -D  <  1 , 1 )  )  /  1C 
9  PENUP 

6222  MG VE  Ll,-(.16*£) 

8223  LABEL  INI  (LI) 

6224  NEXT  LI 

6226  LDIR  0  is  PENUP 

6227  MOVE  D  (1 , 1)  +  .2*  (D  (  300 , 1)  -0  ( 1 , 1 )),-(.  7*E  ) 
o22b  LABEL  "FREQUENCY  (HZ)" 

6230  PENUP  tS  PEN  1 
62  32  LDIR  0 

8234  FOR  L1*-(.6*E)  TO  E+.2*E  STEP  E/5  is  PENUP 
82  3  5  MOVE  D  ( 1 , 1)  - .  09 *  (D  (  20U  ,  1 )  -0  (  1 , 1 )  )  ,L  L 
6236  LABEL  INT(LI) 
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6227  NEXT  LI 
6240  LDIR  0  ,SIN ( 90 ) 

6  241  MOVE  D(1,1)-.1*(D  (300, L)-D  (1,1)  )  1*E 

6243  LABEL  H$ 

6250  RETURN 

6300  REN  **PLOT  IN** 

6201  G CL EAR  §  CLEAR 
6  30  3  H6*D  (1,1).-.  15*  ID  (300,1) -D  (1,1)  ) 

6  304  HS*D(300,1)  +  .1*(D(300,1)-D  (1,1)  ) 

6  305  L*MAX  (AES  (£5  )  ,AtS  (E  3)  ) 

6306  PLOTTER  IS  70S  y  PEN  1 
6312  SCALE  H8  ,K9  ,-L-.  5*L  ,L+  .  25*L 

6212  XAXIS  0,(D(  300  ,1)-C(1,1)  )/lu,C(  1  ,1)  ,D(  200,1) 

6314  YAXIS  C( 1,L)  ,L/5  ,-L-.2*L ,L+ .2*L 

6315  PENUP 

8316  FOR  1*1  TO  200 
6225  PLOT  D  (I  ,1)  ,C( 1 ,  2) 

6330  NEXT  I 

8  3  21  PENUP  y  DEG  §  LDIR  C,SIN(9C) 

6333  FOR  Ll=C  (1,1)  TO  0(  200  ,1)  STEP  (D  (  20  0  ,  1  ) -D  ( 1 ,  1)  )  /  1G 
y  PENUP 

8335  MOVE  LI ,-L- . 2*L 
6336  LABEL  IN.T(Ll) 

6237  NEXT  LI 

8340  LDIR  0  (e  PENUP 

8  341  MOVE  D  (1,1)  +  .2*  (D  (  300 ,1)-D  ( 1 , 1)  )  ,-L-.  4* L 
6342  LABEL  "FREQUENCY  (HZ)" 

6  3^5  PENUP  (>;  PEN  1  e  LDIR  0 

6  352  FOR  LL=-L-.2*L  TC  L+  .2*L  STEP  L/ 5  t>  PENUP 
6355  MOVE  D(1,1)-.C9*  (D(  20G,1)-D  (1,1))  ,L1 
6256  LABEL  INT(Ll) 

3357  NEXT  LI 

8360  PENUP  y  PEN  1  @  LDIR  C,SIN(90) 

6  361  MOVE  D (1 ,1) -. 1* (D (  300  ,l)-0 ( 1,1) )  ,-  (L* . 5) 

6366  LABEL  G$ 

o  360  RETURN 

8400  REM  **CIRCLES* * 

640  1  GCL EAR  €  CLEAR 

8402  L1*MAX(ABS  (E2)  ,ABS  (E5)  ) 

640  2  L*MAX  (AES(E)  ,L1) 

6404  PLOTTER  IS  705  y  PEN  1 

8406  SCALE  1. 25* (-L- . 1*L)  ,  1 . 25* (L+ .  1  *L)  ,-L- . 2*L  ,L+ . 15*L 

6407  XAXIS  0  ,L/5  ,-L-.l*L  ,L+ .  1*L 

64C8  YAXIS  0  ,L/5  ,-L-.l*L  ,L+.  1*L 

6409  PENUP  y  PEN  1 

6411  FOR  1*1  TC  200 

o 4 12  PLOT  C  (I  ,2)  ,L( I  ,  2) 

o  4  1 3  NEXT  I 

8414  PENUP 

6415  LDIF  G  y  PEN  1  y  PENUP 
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6416  MOVE  -<L*.S)  ,-L-.l*L 
6418  LA EEL  H$ 

6420  PENUP  @  PEN  1  @  DEG 

6425  LCIR  C , SIN ( 9Q) 

6426  FCR  Ll— L-.1*L  TC  L+.1*L  STEP  L/  5  §  PENUP 

6427  MOVE  Ll,-L 

6428  LABEL  IN.T<L1) 

6429  NEXT  Ll 

6435  LDIR  0  @  PEN  1  @  PENUP 

8436  FOR  Ll—  L-.1*L  TO  L+.1*L  STEP  L/ 5  @  PENUP 
6437  MOVE  - (.666*L) ,L1 

6436  LAEEL  INT(Ll) 

6439  NEXT  Ll 

8445  LDIR  0,SIN(9G) 

0446  WOVE  - <.95*L)  3*L) 

6446  LABEL  G8 
6450  RETURN 
8549  REN  **NORM** 

6550  N8  fN9*0 

o551  FOR  1*1  TC  10 

6557  WAIT  5u0 

8556  OUTPUT  709  ;"AC3" 

8555  ENTER  7C5  ;  N(I) 

6560  N8=N8+N(L) 

8565  OUTPUT  7G9  ;"AC4" 

0566  ENTER  709  ;  N( I) 

6567  NS*K9+M(L) 

0568  NEXT  I 
8565  N8=N8/ ( I  —  1 ) 

6570  NS— (KS/(I-1)) 

8572  PRINT  USING  6573 
6  57  3  INAGE  3/ 

8575  RETURN 

8600  REN  ** VALUES** 

6601  CLEAR  s  BEEP 

6602  LISP  "I  Ah  COLLECTING  THE  DAIA  NGm" 

6605  IF  L2*l  AND  N  =  2  THEN  6620 

66  06  IF  L2*l  AND  >i*i  THEN  o6  21 
6607  IF  L2*2  AND  K*1  THEN  8622 
o60o  IF  L2*2  AND  N*2  THEN  6623 
6620  F6*F1  s  y*Gl  fe  GCTC  Oo25 

0621  Fo*F0  G  o=w0  £  GOTO  6625 

6622  Fo*FS  s  G=l*6  y  GOTO  0625 

0623  F6*FS  s  W*u9 

8625  FCR  1*1  TC  50 

o628  IF  F  6-4  2*F6/w  +  ( Z  2-2 )  *  F6  /  ( 5  U*G  )  >2  00  THEN  oo31 
do 29  D  (I  ,1)*20Q+(F6-2*F6/v-2ou)  *1/50 
6630  GOTO  0632 

66  31  D  (I  ,  1 ) *F6-22*F 6/G+I* ( 2  2-2 ) * F6/ { 50*v) 

6o32  DS*vALo  (0(1,1)) 
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t>6  3  3  OUTPUT  7 17  ;A$,D$,B$ 

8634  WAIT  150 

8635  OUTPUT  709  ;"AC3VT3" 

8636  ENTER  709  ;  Hti 
6637  D(I  ,2) =H6*S8/N8 
8636  OUTPUT  709  ;"AC4VT3" 

8639  ENTER  709  ;  H9 

6640  D(I  ,3)=H9*£6/N9 

6641  NEXT  I 

66  50  FOR  1=51  TO  2  50 

66  51  D  ( I  , 1) =F6-2*  F  6/0  +  (1-50)  *  2*  F  6/  (  100*0) 

8652  D$*OAL$ ( D { I  ,1) ) 

8653  OUTPUT  717  ;A $,£$,£:? 

0654  WAIT  250 

8655  OUTPUT  709  ;"AC2VT3" 

6656  ENTER  709  ;  HS 
8657  D (1 ,2) =H6  *S8/N8 
6656  OUTPUT  7C9  ;"AC4VT3" 

6659  ENTER  709  ;  K9 

6660  D  ( I  ,  3 ) =H9*S6/N9 
66  o 1  NEXT  I 

66  70  FOR  1=251  TO  300 

06  71  D (I  ,  1) =F6+2*F6/0+ (1-250) * (2 2-2 ) *F 6/ ( S0*U ) 

6672  OS*vAL$(C( I  ,1) ) 
o673  OUTPUT  717  ;A$,D$,E$ 

8674  aAIT  150 

66  7  5  OUTPUT  709  ;  "AC  3VT3" 

6676  ENTER  70S  ;  H8 
bo  77  C (I ,2) =H6*SS/N8 
0676  OUTPUT  709  ;"AC4VT3" 

8679  ENTER  70S  ;  HS 
o680  0  ( I  ,  3 ) =H 9*S6/ N9 
8681  NEXT  I 

£6b2  CLEAR  !g  LISP  "I  AN  FINDING  FIN/NAX  VALUES  &  ASSOC. 
FREwS" 

86b3  GCSUb  870C 


8 684  IF  L2=2  THEN  8687 

8685  PRINT  "GiVAX*"  ;E  ,  FGNX=  " ;  E 1 .  EGRX="  ; E2 
,"FBFX="  ;£<,,  ,"bM="  ;£5,  f BRI  =" ;  L6 

bobb  GOTO  fa  6 97 

86 87  PRINT  "Ri\AX*";E,  c  Ri*.X=  "  ;  £  1  ,  XF>:X  =  "  ;  E2 
FXMX*";E4,  , "  Xi\  I  =  "  ;  E  5  ,  FXNI.*"  ;  £6 
8o  9  7  DISP  "IF  ALL  IS  wELL ,  ENTER  'l'.  ENTER  ' 
E  DATA." 


i 


"  x  i .  V  s  "  •  c  t 

ii  ■>  i 

'2'  TO  RETAK 


8698  INPUT  2 
ab9 9  RETURN 


o7U0  REF.  * *FA Xf-'IN  *  * 

6703  h  ,H  1  ,H2= 100 

6  705  £*D<1U0,2;  (e  El  ,£4  ,E6=D  (10o  ,1)  «  E  2  ,  E  3  ,£  5=D  ( 10  8 , 3) 
J/Gb  FOR  I=1G1  TO  20o 
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8710  IF  E(I,2)<D|H,2)  THEN  6725 
6711  IF  C(I,2)*C(H,2)  THEN  8720 
8720  E*D  (1,2)  fc  E1*D(I,1)  t§  E2»C(I,3)  £  H*I 

6725  IF  C  (I  ,3)  <D  (Hi  ,3)  THEN  8740 

6726  IF  C(I  ,3)*D  (Hi, 3)  THEN  67  35 
6735  E3*D(I,3)  @  E4-0(I,1)  6  H  1*1 
6740  IF  C ( 1 , 3 ) >*C (H 2 , 3 )  THEN  6751 
6745  H2*I  @  E5*D  (1,3)  e  £6*D(I,1) 

6751  NEXT  I 

6760  J  1  *D  ( H  ,  2 )  U2=(D  <H+l,2)-0  <H-1  ,2)  )/2  f  J3=(0(H  +  i,2) 

+0(H-l,2).-2*dl)/2 

8762  F6*-(J2/(2*J3) ) 

6  7o  4  A*J  1+J2+F6+J  3*F6“  2 
o  770  Fb*C  (H  ,1)+F6*  (D(H  ,1)  -D(H-1,1)  ) 
b  7  7  5  E1*F6  tf  E*A 
6766  D=ABS (E  3-E  5 ) 

6765  IF  M*2  THEN  87  5  3 
6750  IF  L2= 1  THEN  6752 
b/yl  C 3*C  @  GCTC  b 7 56 
8792  D1=C  e  GOTO  6756 
o  75 3  If  L2= 1  THEN  8  79  5 
6754  D4*C  £  GOTO  6756 
6795  C2=C 

o7So  PRINT  "  D 1*"  ;  D 1 ,  ,  "C2*  C2  ,  ,"D3*";D3  ,  ,"0  4*"  ;  04 
o 7  5  5  RETURN 
9  oCO  REP.  **CALC-A*  * 

5001  PRINTER  IS  701,76 

5002  PRINT  uj  PRINT 

9003  PRINT  "VALUES  MEASURED  IN  AIR  FPCM  M;8S;"  DATA" 

50u4  PRINT 

5011  PRINT  “ BLOCKED  CAPACITANCE*  H;A£S(Co);"  FARADS  AI  " 
;F ; "HZ" 

9bl2  PRINT  "ELECTRICAL  V.UALIT*  FACTCR  =  "  ;w0 
y 01 3  PRINT  "MECHANICAL  vOALIIY  FACTCR  =  " ; ,o 

5014  PRINT  "ELECTRICAL  RESONANCE  *";r0;"  Hi." 

5015  PRINi  "MECHANICAL  RESONANCE  =  ";cc;”  LZ  " 

5  0  i  o  PRINI  "DYNAMIC  ELECT  rCMECHAN  I  CAL  COUPLING  COEFFICIE 
NT  *  "  ;  K 

3017  PRINT  "STATIC  CCUPLING  CCtf  r  ICIci/i  *  "  ;  *  1 
502  1  PRINT  "BLOCKED  RESISTANCE*"; Po;  “CrN-S" 

5031  PRINT  uSING  5032 

5032  IPAGE  2/ 

9 0 5u  RETURN 

5  lo  0  REM  **CALC-N*  * 

5101  PRINTER  IS  701,76 
9  Lo  2  PRINT  ^  PRINT 

5  lo  :  PRINT  "VALUES  t  CR  i>  E  A  S  U  R£i\  L  NT  £  IN  ,vAT£R  FRCR  ";£8; 
"  DA iA" 

5 lu  4  PRINT 

5lU  PRINT  "  BLOCKED  CAPACI LANCE  =  ";AL3(C0);"  tARAOS  AT 
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" ;  F  ;  "HZ" 

a  112  PFINT  "ELECTRICAL  wUALITY  FACTOR  *  M;CJ 
9112  PRINT  "MECHANICAL  QUALITY  FACTOR  *  ";v£ 

S 114  PRINT  "ELECTRICAL  RESONANCE  *  ";F1;H  HZ" 

9115  PRINT  "MECHANICAL  RESONANCE  =  " ;F9;"  HZ" 

9116  PRINT  "DYNAMIC  ELECTRCMECHAN ICA.  COUPLINC  COEFFICIE 
NT  *  " ; K 

9120  PRINT  "MECHANICAL  POWER  UTILIZATION  FACTOR" ; AES ( P6 ) 

9121  PRINT  “  STATIC  COUPLING  COEFFICIENT*" ;Kl 

9122  PRINT  "ELOCKED  RESISTANCE  =" ; R9; "OhM-S " 

9122  PRINT  USING  9133 

912  3  I  PAGE  3/ 

9150  RETURN 
S 200  REM  **CALC-A/M** 

9201  PRINTER  IS  701,76 

9202  PRINT  ig  PRINT 

92 02  PRINT  "VALUES  MEASURED  IN  AIR  FROM  ";S9;"  DATA" 

9204  PRINT 

9211  PRINT  "BLOCKED  INDUCTANCE*  " ; ABS (LC )  ; "  HENRIES  AT  " 
;F6 ; "HZ" 

9212  PRINT  "ELECTRICAL  wUALITY  FACTOR  =  "  ;  *0 

9213  PRINT  "MECHANICAL  wUALITY  FACTOR  =  " 

9214  PPINI  "ELECTRICAL  RESONANCE-  =";Fu;"  HZ" 

9215  PRINT  “MECHANICAL  FtSCNANCE  *  " ;Fb;"  HZ" 

92 lo  print  "DYNAMIC  ELECTROMECHANICAL  CCLPLINC  COEFFICIE 
NT  *  "  ;  K 

9217  PFINT  "ELCCKEC  RESISTANCE*  “ ;F;"CHNS" 

9231  PRINT  USING  9232 
92  22  IMAGE  2/ 

9250  RETURN 
9 20 o  FEM  **CALC-N/M** 

93E1  PRINTER  IS  701,76 
9  3C2  PRI^l  ip  PRINT 

9  30  3  PRINT  "VALLES  FCF.  MEASUFENERTS  IN  a  AT  ER  FPCM  "  ;5$; 
"  DATA" 

9  204  PRINT 

£311  PRINT  " ELCC RED  INDUCTANCE  *  ”;AcS(Lo);"  HENRIES  AT 
" ;Fo;  "HZ" 

a  312  PRINT  "ELECTRICAL  wuALITY  FACTOR  =  ”;wl 
9312  PRINT  "MECHANICAL  wUALITY  FACTOR  =  " ; w £ 

9214  PRINT  "ELECTRICAL  RESONANCE  =  "  ;  F 1 ;  "  HZ" 

9215  PRINT  "MECHANICAL  RESONANCE  =  ";F9;"  HZ" 

S  3  16  FhlrtT  "DYNAMIC  ELECTRCMECiiAN  I  CAL  COUPLING  COEFFICIE 
NT  *  " ; K 

9  2 lo  PRINT  "POTENTIAL  EFFICIENCY  =  ”;P5 

9  219  PRINT  "FFE0ULNCY  OF  OPTIMUM  EFFICIENCY  =  "  ;  o ;  "  HZ" 
9  320  PRINT  "MECHANICAL  PC  ,v£F  UTILIZATION  FACTOR  "  ;AcS  (Po  ) 
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S  322  PFItil  USING  5323 
s*  323  IN£GE  3/ 

9  3  SO  PETGEN 

5556  DISP  "  *IK£  END" 

5955  END 


no 


1  REM  **PICK£W** 

2  OPTION  EASE  1 

3  SHORT  C0,Gl,G8iQ9»Q,E,E2,£3iE5,R8,R7,X2,X3,X5,X6,X7,X8 
,G1 , G8  ,  E2  ,  B  3  ,  £5  ,  £6  ,£7  ,£8,K,K1,P5,P6,W 

7  SHORT  ArF2,F3,F4tF5tH2,H4,Jl,J2,J3,J4,G,E,RtX 

8  INTEGER  M,  Y  ,  Y I ,  I  ,Y2  ,J  ,  22  ,  Z  3  ,Hi  ,H2  ,L  2  ,S  1 

10  DIM  B$  {  2)  »A  $  ( 2  ]  »>T$  f  1]  ,D  ( 300 , 3)  ,A  1$  [4]  ,D$  [  30C  J  ,L$.  [20 )  , 
RRfll  ,I$[1J  ,S$ [ 10]  ,C$  ( 5]  ,G$  [2  2]  ,H$.  (2  3] 

11  REAL  S  ,N  ,N8  ,N9  ,R0  ,R9  ,X0  ,X9  ,L  ,C0  ,LC  ,0  ,P1  ,H8  ,H9  ,P,W0,G0 
iG9  ,E0  ,89  ,N  ( 10)  ,M ( 1 Q. )  ,F  0  ,F1  ,F  6  ,  F6  ,F  9 

14  REAL  S8  ,E4,E6  ,L5  ,D  ,Dl,D2  ,D2  ,D4 

15  X6,Dl,D2,D3,D4  =  0 

16  C,E2,E3,E5,R8,R7,X2,X3,X5=0 

17  Z3,A,F2,F3,F4,F5,B6,£7,£8=C 

18  S8  ,E4,E6  ,M,Y,Yl,I  ,Y2,d  ,Z2=0 

19  K  ,Kl  ,P5  ,P6  ,W  ,C0  ,Q1  ,£8  ,C9  ,£=0 

20  CO  ,L0  ,0  ,P1  ,B  ,  H6  ,H9  ,P  ,W0  ,S  =0 

21  N  ,N8  ,N9  ,R0  ,R9  ,X0  ,X9  ,G0  ,G9=0 

22  BO  ,B9  ,Hl  ,H2  ,F0  ,F  1  ,F6  ,F8  ,FS*0 

24  DISP  "THIS  PROGRAM  IS  DESIGNED  FOR  AN  EXPERIENCED  CPE 
RATCR.  THE  BANDWIDTH  FOR  DATA  COLLECT . " 

25  DISP  "AND  FREC.  TO  GET  BLOCKED  DATA  *ILL  BE  THE  ChCIC 
E  OF  THE  OPERATOR.  "HYDRA 2'  PROGRAM  IS" 

26  DISP  "DESIGNED  TC  COVER  1C  X  6W  AUTOMATICALLY.  HIT  "C 
CNT '  IF  READY  TC  FRCCEED." 

27  PAUSE 

28  CLEAR  @  BEEP 

.35  DISP  "ENTER  TYPE  OF  TRANSDUCER  TC  BE  MEASURED?  ("m"  FC 
R  MAGNETIC  COUPLING  CR  "E"  FOR  ELEC.)" 

42  INPUT  T$@  D ISP  "TC  GET  A  COMPLETE  SET  OF  DATA  YOU  NEE 
D  MEASUREMENTS  IN  BOTH  AIR  AND  aATER.  ” 

45  DISP  "DC  AIR  FIRST  . IN  aHAI  MEDIUM  APE  YOU  OPERATING?  ( 
ENTER  "!'  FOR  AIR  OR  "2"  FOR  WATER)” 

46  INPUT  Kg  CLEAR  BEEP 

47  IF  Ml =2  THEN  49 

4o  CS="AIR"  0  GOTO  51 
49  C$= "WATER" 

51  PRINTER  IS  2  ^  PLOTTER  IS  1 

55  DISP  "SET  THE  CRANE! Z  FOR  ACMI TTANCE  .SET  ON  LOWEST  SC 
ALE  THAT  WON'T  PEAK  DURING  RUN.  iCCNT ) " 

6C  PAUSE 

6  5  PRINT  («  PRINT  @  PRINT 
70  PRINT  "ADMITTANCE  IN  ";C$ 

75  GCSUE  7000 

60  IF  M*2  THEN  SO 

85  F0-F6  <e  00=0  is  GOTO  100 

90  F 1=  F6  ^  s.l*V 

100  GCLEAR  @  CLEAR  BEEP 

1 C 5  DISP  "SET  DRANETZ  FOR  IMPECANCE-(SCALE  NEEDED ) -HIT  C 
CNT  ” 


110  PAUSE 

115  PRINT  @  PRINT  @  PRINT 
120  PRINT  "IMPEDANCE  IN  ";C$ 

125  GOSUB  7000 

130  IF  M=2  THEN  140 

135  F8*F6  e  Q8*G  @  GOTO  141 

140  F9*F6  @  u9=G 

141  CLEAR  £  BEEP 

142  PRINTER  IS  701,76  &  PRINT  USING  144 

144  IMAGE  2/ 

145  PRINT  "F  (  Y-AIR )  *  "  ;  INT  (F  0)  ,  ,"F  (fc-WAT)  *  "  ;IUT  (F  3)  ,  F 
( Z  -  A I R )  =  "  ;  IN.T  ( F8)  ,  F  (Z.-WAT)  ;  INT'(F&) 

146  PRINT  "C(Y-AIR)=  “  ;  GO ,  w  ( Y-NATEF)  “  "  ;G1 »  »  HW‘  (  Z-AI R )  * 
";u8  ,  ,"C(Z-WATER)  =  ";w9 

147  PRINTER  IS  2 

148  CLEAR  @  DISP  "THESE  ARE  THE  INITIAL  ESTIMATES.  (HIT  ' ' 
CCNT'  WHEN  READY  TO  PROCEED) “ 

149  PAUSE 

150  GOSUB  2501 

152  CLEAR  6  BEBP 

153  DISP  "TO  COLLECT  ADMITTANCE  DATA  (FOR  ELEC.  COUPLING 
)  ENTER''!'.  ENTER  "2"  FOR  IMPEDANCE" 

154  INPUT  L2@  PRINTER  IS  2 

155  IF  L2=2  THEN  158 

156  R$="G"  I  $  - "  B "  t?  S  $="  ADMITTANCE"  (‘  H$  =  "CONDUCTANCE 
(MI CRCMHCS) "  £  GS-"SUSCEPTANC£  (MICRCMHCS ) " 

157  GCTC  159 

158  R$=  "R "  C=  I$="X"  t  S$*"  IMPEDANCE"  e  G$="  REACTANCE  (OH 
MS)"  &  H$="RESISTANCE  (OHMS)" 

159  IF  L2= 1  THEN  176 

16C  CLEAR  @  DISP  "SET  ON  Z  AND  ENTER  SCALE  FACTOR  TO  COL 
LECT  DATA . " 

161  INPUT* S8@  PRINT  S$;"  IN  ";C6  b  PRINT  "SCALE 
OHMS" 

162  CLEAR 

164  GOSUB  8600 

165  IF  Y  =  2  THEN  152 
loo  IF  M*2  THEN  169 

167  R8=£  (§  X3  =  E2  <s  X6=E  3  <a  F4=E4 
16a  X6=E5  «  F5=E6  (g  F8  =  El  fe  GOTO  19b 

169  F7*E  6  X2=t 2  v  X7=E3  t  F4=E4 

170  X5*E5  e  F  5=£6  *  F5-E1  &  GOTO  19b 

176  CLEAR  b  BEEP 

177  DISP  "SET  ON  Y  &  ENTER  SCALE  FACTOR  IN  MICRCMHCS  TO 
COLLECT  DATA." 

178  INPUT  S6@  PRINT  S$;"  IN  "  ;C$  C°  PRINT  "SCALE  =";S6;" 
MICROMHOS" 

179  GOSUB  8600 

160  IF  Y»2  THEN  152 

161  IF  M-2  THEN  185 
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Ie2  G6*£  8  8 3*E 2  8  86«E3  8  F2«£4 
183  B6*E5  8  F3-E6  8  FG»E1  8  GCTC  196 

165  G7»E  8  B2-E2  8  67-E3 

166  F2«£4  8  B5«E5  8  F2*E6  8  F1*E1 
1S6  PRINT  USING  199 

199  I  PAGE  2/ 

200  BEEP  8  DISP  "DC  YOU  DESIRE  A  LIST  OF  DATA  NEAR  RESON 
ANCE?  (10C  PO INIS)  (1»YES  ,2«NC) " 

201  INPUT  Y 

202  IF  Y»2  THEN  300 
204  IF  L2*2  THEN  240 

210  PRINT  "ADMITTANCE  DATA  FCR  ";C$;"  IN  MICF.OMHOS"  8  PR 
INT  USING  215 

215  IMAGE  2/  FREQUENCE" , 4X  , "REAL "  ,6X  , "IMAGINARY" 

216  GOTO  265 

240  PRINT  "IMPEDANCE  DATA  FOR  IN  OHMS"  8  PRINT  US 

ING  215 

265  C-CSUB  8100 

3GU  CLEAR  8  BEEP  8  DISP  “DO  YOU  DESIRE  A  PLOT  OF  DATA?  ( 1 
=  YES  ,2*NC) " 

30  5  INPUT  Y@  GCLEAR  8  CLEAR 
315  IF  Y=2  THEN  450 

320  DISP  "ENTER:  1*  G/R  MS.  FREw;  2*  E/X  MS.  FREG;  3“  L 
MS.  G/X  MS.  R;  4=  END  PLOTTING  LOOP." 

325  INPUT  Y 
320  IF  Y»2  THEN  365 
335  IF  Y»3  THEN  425 
240  IF  Y=4  THEN  450 
345  GCSUE  6200 

356  LCIR  0  8  PEN  1  8  PENUP 

357  MOVE  0(1,1)+. 2*  (D  (300,1) -0  (1,1))  ,-(E* . 6) 

360  LABEL  R8;"  MS.  FREQUENCY  IN  " ;C$ 

361  GOTO  320 
365  GOSUB  8300 

397  LCIR  0  8  PEN  1  8  PENUP 

298  MOVE  D(l,l)  +  .2MD(3Go,l)-0(l,l))  ,-L-.5*L 

400  LABEL  1$;"  VS.  FREQUENCY  IN  ”;C8 

401  GCTC  220 

425  GCLEAR  8  CLEAR  8  GOSUE  64vC 
421  FEN  1 

442  LCIR  0  8  PENUP 

443  MOVE  -(. 5*L) ,-L-.2*L 

445  LABEL  "INPUT  ELECTRICAL  ";S6;”  PLOT  (";I $  ;  " vS. " ; E$; " 
)  IN  " ; CS 

446  GOTO  320 

450  IF  L2*2  THEN  645 

451  IF  T$*"E"  THEN  525 

452  GCLEAR  8  CLEAR  8  DISP  "I  AM  CALCULATING  MCTICNAL  DAT 
A" 

465  CC-E*.U00001/(2*PI*F6) 
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466  L0*1/ ( (2*PI*F6)~2*C0) 

47C  FOR  1*1  TO  300 

475  D(I  ,2)*D(I  ,2).-G 

480  D  (I  , 3) *D  (I  ,  3).-D  ( I  r  1 )  *B/F6 

485  NEXT  I 

490  GOSUE  8700 

500  GOTO  542 

525  GCLEAR  @  CLEAR  §  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 

A" 

532  C0*B*.GQ0001/( 2*PI *F ) 

535  FOR  1  =  1  TO  300 

536  D ( I ,2 ) *D  ( I  ,2).-G 

537  D  ( I  #  3 )  *D  1 1  f  3  ).-C  ( I  #  1 )  *B/F 

538  NEXT  I 

540  GOSUE  8700 

542  BEEP  <§  DISP  "DO  YOU  DESIRE  A  LIST  OF  MOTIONAL  DATA  N 
EAR  RESONANCE? ( 1=YES  ,  2=NC)" 

54  3  INPUT  Y 

545  IF  Y=2  THEN  604 

575  PRINT  "MOTIONAL  ADMITTANCE  DATA  IN  " ;  C$  §  PRINT  OSIN 
G  £80 

560  IMAGE  2/  /’FREQUENCY"  ,4X,"REAL"  ,6X  /"IMAGINARY" 

555  GCLEAR  g  CLEAR 

6 CO  GOSUE  8100 

604  CLEAR  @  BEEP  @  BEEP 

6C5  DISP  “KANT  A  PLOT?  (  1*YES  ,2=*NO)  " 

610  INPUT  Y2 

615  IF  Y 2=2  THEN  755 

616  PLOTTER  IS  705  @  PEN  1 
620  GGSUB  8400 

628  LDIR  0  S  PENUP  @  PEN  1 

629  MOVE  -  (. 5*L)  ,-L-.2 *L 

630  LABEL  "MOTIONAL  ";S$;"  PLOT  FOR  " ;C$ 

642  GOTO  755 

645  CLEAR  (?  GCLEAR  £  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 

A" 

646  IF  T$*"K"  THEN  65o 

647  FOR  1*1  TO  300 
b46  D  ( I  ,2)  *C(  I  ,2)  -F 

649  C0*1/ (2*PI*F*X) 

650  LG* 1/ ( ( 2*  PI  *  F6 ) “ 2*C  0) 

651  D  ( I  ,  3)  *0  ( I  ,  3  ).-2*P I *D  ( I  ,  1)  * L0 

652  NEXT  I 

653  GO SUB  8700 

654  C0*A£S (CO )  @  L0*ABS (LO) 

655  BEEP  £  DISP  "DC  YOU  KANT  A  LIST  OF  THE  DATA  NEAR  RES 
CNANCE? ( 1*YES ,2*NC) " 

656  INPUT  Y 

657  GOTO  670 

658  L0*X/(2*PI*F6) 


662  FCR  I«1  TO  300 

663  D ( I #2 ) *D ( I  ,2).-R 

664  D  ( I  i  3 )  *0  (1,3 ).-D  (1 ,1)  *  2*P  I  *  LO 

665  NEXT  I 

666  GOSUB  8700 

667  BEEP  @  OISP  "  DO  YOU  WANT  A  LIST  OF  DATA  NEAR  RESGNA 
NCE?(l*YES,2«NCt)" 

668  INPUT  Y 

670  IF  i»2  THEN  706 

685  PRINT  "MOTIONAL  IMPEDANCE  DATA  IN  ";C$  <s  PRINT  USING 
666 

6fa6  IMAGE  2/  , "  FREQUENCY"  ,  4X  , "REAL "  ,bX  ,  "  IPAGINARY " 

700  GCLEAR  CLEAR 

705  GOSUB  8100 

706  CLEAR  @  BEEP  BEEP 

710  DISP  "WANT  A  PLOT? ( 1=YES  ,2-NC) " 

715  INPUT  Y2 

720  IF  Y2=2  THEN  755 

721  PLOTTER  IS  705  @  PEN  1 
725  GOSUB  84G0 

737  LDIB  C  @  PENUP 

738  MOVE  -  (. 5*L)  ,-L-.2*L 

740  LABEL  "MOTIONAL  PLOT  IN  " ; C$ 

755  GCLEAR  §  CLEAR 

756  PRINT  USING  757 

757  IMAGE  3/ 

SCO  DISP  "I  AM  DOING  CALCULATIONS  FOR  YOU.  PLEASE  BE  ffiT 
I  ENT." 

603  IF  M»2  THEN  1040 
610  IF  TS*="M"  THEN  1000 
625  IF  L2*2  THEN  660 
640  00»E1/A£S (E4-E6) 

045  C1=ABS(E2/(2*PI*FQ) ) 

655  GOTO  865 

06  0  g«»E 1/A BS (E4-E6) 

665  Cl*AES( 1/ (2*PI*FS*E2) ) 

oo5  K* 1- (F0/F6) *  2 

855  K1C  1/(ABS(C0)+C1) 

5:0  GOSUB  9000 
931  PRINTER  IS  2 
540  GOTO  1350 
1CU0  K*l- (F6/F0) “2 
1005  IF  L2«2  THEN  LU25 
1010  wO*£l/ABS(E4-Eb) 

1015  GOTO  10  30 

1025  w6«El/AES (E4-E6) 

It 30  GOSUB  5200 

10  31  PRINTER  IS  2  is  GOTO  135u 

1040  IF  X6*"M"  THEN  1150 

lt4 1  BEEP  IS  DISP  "IF  YOU  HAvE  DATA  FCR  AIR,  ENTER  1;  ELS 
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E  2" 

1042  INPUT  Y 

1043  RAD 

1044  IF  L2»l  THEN  1110 

1045  IF  Y»2  THEN  1062 

1060  P6*D4* (D3-D4) /  (C3*R7) 

1061  0*  F  1 

1062  Q 1*E1/ABS (E4-E6)  @  K=  l-( F  1/F 9)  “2 
1070  GOSUE  9100 

1076  PRINTER  IS  2  g  GCTC  1350 
1110  IF  Y*2  THEN  1135 
1126  foG*F] 

1130  P6*D2* (D1-D2)/ (D1*G7) 

1135  C1*£1/ABS(E4-E6) 

1140  K*l- (F1/F9) "2 

1145  GCSUE  9100 

1146  PRINTER  IS  2  fc  GOTO  135u 

1150  RAD 

1151  DISP  “IF  YCO  HAVE  DATA  FOR  AIR,  ENTER  1;  ELSE  2" 

1152  INPUT  Y 

1153  IF  L2*l  THEN  1160 

1154  IF  Y*2  THEN  1162 

1155  B=.5*ACS( (R7-R9)/D4) 

1156  H6*SIN<B)*2  $  H9*C0S(0)“2 

115  7  P5*(  ( 1+D 3/R0 *H  9 )  “.5-(  1-D3/P0*H6)  “ .  5 )/  (  (  1+D 3/R0 *ri 9)  * 
.  5+  ( 1-D  3/H0*Hd)  **  .5} 

1159  P=Q4*SIN  (2*E)  /  (4*R9*<*9)  t  bG* F 9* ( 2* P + ( r '2+ l ) “  . 5 ) 

1160  P6=D<»* (D3-D4) / (D  3*R  7) 

1162  C9*E  1/ABS  (E4-E6)  <?  K«l-(r  VFl)‘l 

116  3  PRINT  "B»"  ;  E , "  h*  " ;  K  K  i»"  ;  K  1 ,  “  he-"  ;  no  ."  n  9=  "  ;  n9 
1164  GCSUB  9300 

llbd  PRINTER  IS  2  «  GOTO  1350 
1160  RAD 

1182  IF  Y*2  THEN  1210 
1164  E*.5*ACS( (G7-G9)/C2) 

1166  H8*SIN(S)~2  g  H5»CCS(b)*2 

1166  P  5*  (  ( 1+C1/GG*H9)  5- ( 1-D  1/Go* Ho  )  '  .  5  )  /  (  (  i+C  1/G0 *U  5  ) 

.  5+  ( 1-C  l/CO*Hd)  “  .5) 

1190  P.*G2*SIN  (2*H)/ (4*G5*wl)  b  aC*  F 1  *  (  2*P  +  ( P  "2  1 )  ”  .  5  ) 
1195  P6*D2*  (C1-D2)  /  (C1*G/ ) 

1210  w9*El/A£S  (E4-E6)  a  R«L-(F  9/F  1)  “2 

1211  PRINT  "£*"  ;xs  ,“K*  H ;  ft  a  1-  M ;  K  1  ."K8-"  ;uo  ,"h9-";t:9 
1215  GCSUB  930o 

1225  PRINTER  IS  2  J  GOTO  1350 

1350  DISP  "DO  YCU  DESIRE  ANOTHER  IN  THIS  >.ECIuR>(TO  GET 
THE  OTHER  TYPE  DATA?)  ( 1*YES  ,2*f.G)  " 

1355  INPUT  Z1 
1357  PLOTTER  IS  1 

1360  IF  21*1  THEN  152 

1361  PRINT  US  INC  136  2 
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1362  IMAGE  3/ 

1365  CISP  "DO  YOU  DESIRE  A  RUN  IN  ANOTHER  MEDIUM?  ( 1*YES  , 
2*NO)  " 

13/0  INPUT  21 
1375  IF  21*1  THEN  45 

1360  DISP  "JUST  FOR  THE  RECORD,  INPUT  WATER  TEMP,  AIR  T £ 
MP,  TRANSDUCER  SER  NO,  MODEL  NO.," 

1361  DISP  "AND  TYPE  (E  OR  M)" 

1365  DISP  "INPUT  "O'  IF  INFO  IS  UNKNOWN" 

1390  INPUT  W  , A  ,S  ,N  ,L$ 

1351  PRINTER  IS  701,76 

1400  PRINT  "aATER  TEMP*"  }  w  ,  ,"AIP  TEMP*"  ;A  ,,"  SER.  NO.  =  "?S 
,,"  MODEL  NO.*"  ;N  ,  ,"  TYPE  ";L6 
1410  GOTG  9996 

2501  REM  *3 HUNT* 

2502  CLEAR  £  DISP  "INPUT  A  PERCENTAGE  OF  RESONANT  r  REG. 

TO  USE  TO  GET  S HUNT  (BLOCKED)  VALUES." 

2503  DISP  "INPUT  A  VALUE  FROM  0  TO  1." 

2504  INPUT  £  1 

2  513  CLEAR  g  BEEP  DISP  "WE  ARE  FINDING  SHUNT  VALUES  FO 
R  G/B  CR  R/X." 

2515  A  $* "FR "  &  £  $.  *  "  H  Z  "  &  Sl*l 

2519  IF  T  $  * "M "  THEN  2525 

2520  IF  M* 2  THEN  2522 

2521  F=F0*B1  (g  GCTC  2535 

2522  F*F  1  *B 1  e  GOTO  25  35 
25  25  IF  M* 2  THEN  2527 

2526  F*Fa-Bl*F6  fc  GOTO  2535 

2527  F  *F  S-B 1*F9 

2535  D$*VAL$lF.) 

2537  OUTPUT  717  ;A$,D$,C$ 

2536  WAIT  1000 

2540  CISP  "SET  DRANET2  FREw  .  SCALE  TO  CCvEK  FREw  .  ON  S  YN 
ThESIZER .ZERO  METERS. SET  ON  PS.(CGNT)" 

25^  1  PAUSE 

2542  CISP  "HERE  wE  GET  THE  NORMALIZATION  FACTORS." 

2543  OUTPUT  7U»  ;"vll" 

2546  GGSUB  6549 

255U  G,U,R,X*0 

2551  A$*" FR"  *  E$="Hz"  *  D$=vA LS(F) 

2553  OUTPUT  717  ;AS,CS,E$ 

2561  CLEAR  «*  CISP  "SET  DRANETZ  ON  i  ,S£T  SCALE  , SET  CN  NCR 
«  F OF  FILTER  +  PHASE .PLUG  IN  TRANSDUCER . ” 

2562  BEEP  «  DISP  "ENTER  SCALE  FACTOR  IN  i>  ICRCi-.r.CS . SAwE  S 
CALE  AS  FCR  DATA. (w£  GET  SHUNT  G/B)" 

2564  INPUT  E6&  CLEAR  fc  CISR  "I  Am  WORKING  TO  CE'I  SkLM  v 
ALuES"  p  EE&P 

2565  FCR  1*1  TO  10 

256o  GuaPUT  7  L  7  ;A6,D$,£$ 

2567  >.  A 1 1  lOOU 
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25ob  GUTPul  709  ;"AC3" 

2569  ENTER  70S  ;  M(I) 

2570  G*G+fc(I) 

2571  OUTPUT  709  ;“AC4“ 

2572  ENTER  70S  ;  N(I) 

2573  E*G+N(I) 

2574  NEXT  I 

2575  G*G/(I-1)  e  £*E/(I-j) 

2576  G»G*S8/N6  (s  £*B*S6/N9  &  BEEP 

2576  DISP  “SET  CRANETZ  CN  Z;SET  SCALE  rCR  :AX  RESPONSE .b. 
OvE  T  PANS .  INPUT.  ENTER  SCALE  FAC1CP" 

257S  INPUT  Sole  CLEAR  y  DISP  "I  AN  rCPICIKG  TC  GET  SHUNT  V 
ALOES “ 


2560 

EOF  1*1  TC  lu 

25dl 

OUTPUT  717  ;A$,D$,b$ 

2562 

kAIT  1000 

2583 

OUTPUT  709  ;  “AC  3“ 

2564 

ENTER  70S  ;  V.  (I) 

2585 

R*R+M ( I ) 

2566 

OUTPUT  709  ; "AC 4 " 

2  56  7 

ENTER  7uS  ;  N(I) 

256b 

X*N ( I) +X 

256  S 

NEXT  I 

25SC 

R*R/(I-1)  @  X*X/  (1-3) 

25S 1 

R*R*Sb/N8  B  X*X*S &/N S 

2592 

IF  T  THEN  2610 

2595 

IF  N*1  THEN  2b Go 

2596 

RS*R  vr  X9=X  t  C9=C  B 

B  9=B 

$  GOTO  2601 

2  6CC 

GG*G  t  iuU*6  b  Ru=R  t 

X0  =  X 

2601 

PRINT  “SHUNT  VALUES" 

2602 

PRINT  USING  2603 

4.  b  0  j 

I  RAGE  1/ 

2  60  5 

PRINT  "GO*";G0 , f“60=" 

;B 0 , 

R0*  "  ;RO  ,  .“XC*"  ;a0,  ."GS*" 

;GS , 

II  -  0  =  "  .  t  C  '<  r  t-H.nC  <• 

,  Cii—  ,  t,  y  I  I  K5-  |B-|  1 

X9  =  " 

;  xs 

26  06 

PRINT  b  GCTC  2630 

2b  10 

IF  S 1*2  THEN  262  u 

2611 

0  1*G  t  O  <.*E  t  0  3*R  B 

0  4=X 

26 12 

IF  R.*2  THEN  2b  15 

2613  F.-Ffa  +  L  1  *F  &  «•  Sl*2  B  GOTO  2550 
2615  F=FS+B1*F9  g  Sl*2  fc  GOTO  2550 

2620  G*.5*(G+U1)  e  o=.5*(b+J2)  fc  R=.5*(R+«j3)  b  A  =  .5*(J*»  + 

2)  (e  GCTC  2595 

2630  GCLEAR  «?  CLEAR  t  BEEP 

2635  CIS?  "ENTER  "l'  IF  ALL  IS  NELL;  '2'  IF  TOO  NEED  A  R 
EPEAT 

2640  INPUT  Y 

2645  IF  Y*2  THEN  2501 

2650  RETURN 

7000  REN  **F-V** 

7001  GCLEAR  CLEAR 
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7005  DISP  "ENTER  LOWER,  UPPER  FREC.  FOR  SnEEF.IF  RESCNAN 
C£  UNKNOWN  USE  10000,200000.  (F1,F2)” 

7010  INPUT  F4,F2 

7011  LISP  “ENTER  RMS  VOLTAGE  NEEDED.  (30C0  M'v  FOR  LOW  FRE 
L.  AS  LOW  AS  2C0C  MV  FOR  HIGH) " 

7012  INPUT  Afc  A1$*VAL$(A) 

7020  CLEAR  6  BEEP 

7025  A$*"AM"  @  BS*“MR“ 

70  35  OUTPUT  717  ;A$,A1S  ,B$ 

7040  OUTPUT  70S  ;“AC3VT3" 

7045  F3»CEIL( (F2-F4)/30C) 

7060  CLEAR  @  DISP  "I  AM  WORKING  TO  GET  THE  SPECTRUM  LATA 
FOR  YOU, .THEN  WILL  SHOW  &  MAKE  A  PLOT" 

7065  FOR  I»1  TO  300 
7070  D(I,L)»F4+I*F3 
7075  A$*"HZ,“  B$*"FR" 

7063  D$*vAL$ ( C I I ,1) ) 

7065  OUTPUT  717  ;S$,D$,A$ 

7090  OUTPUT  709  ;"AC3VT3" 

7095  WAIT  100 

7100  ENTER  709  ;  C(I,2) 

7105  NEXT  I 

7110  F5*F4- . 1* (F2-F4 )  @  F6=F2  + . 1* (F 2-F 4) 

7115  F6=F2-h  .  1*  (F2-F4) 

7119  PLOTTER  IS  1  (?  G  CL  EAR  £  CLEAR 
7125  SCALE  F5,F6,-.l  ,1.2 

7120  XAXIS  0 , 2000  ,F4  ,F2 
7125  YAXIS  F4,.l,-.1,1.2 
7140  FOR  1*1  TO  200 
71«5  PENUP 

7150  PLOT  0(1,1)  ,0(1 ,2) 

7155  NEXT  I 

7156  GRAPH  e  COPY 

7160  DISP  "NEED  ANOTHER  RUN?  (  1=  YES  ,2  =  NG  )  “ 

7161  INPUT  Y 

7162  IF  Y*1  THEN  7005 

7163  CLEAR  (§  G CLEAR  (s  BEEP 

7165  LISP  “ENTER  DECISION  POINT  FOR  AMPLITUDE  (0  TO  1.2) 
. (LESS  THAN  THE  MAX  DISPLAYED) .“ 

7170  INPUT  S$  CLEAR  BEE P 

7160  PRINT  "AMPLITUDE  IN  vOLTS";"  FREQUENCY" 

7161  PRINT 

7190  FOR  1*1  TO  200 

7195  IF  D ( I  ,2) <S  THEN  7210 

7200  PRINT  USING  72U5  ;  O (I  ,2) ,D  (I  ,1) 

7205  IMAGE  IX  ,C .  DC  CLDC  CL  ,  1CX  ,LDLLlC  .DC 

7210  NEXT  I 

7215  FRINT  USING  7220 

7220  IMAGE  3/ 

7225  CLEAR  *  BEEF 
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7226  DISP  "IF  YOU  ARE  READY  TC  CONTINUE  ENTER  '1';  OTHER 
WISE  ENTER  '2  ' 

7227  INPUT  Y 

7228  IF  Y=2  THEN  7I6C 

72  35  DISP  "WHAT  IS  CENTER  FREQ,  HALFWIDTH  TC  CONSIDER?" 
7236  INPUT  F,N(?  CLEAR  €  BEEP 

7240  DISP  "I  AM  WORKING  TO  GET  GOOD  FREU  AND  0  DATA" 

7241  FOR  1*1  TO  300 

7242  D(I,1)*F.-N+I*N/150 

7243  A$*"HZ"  @  B$="FR" 

7245  D$*VAL$(D<I,1)  ) 

7255  OUTPUT  717  ;E$,D$,A$ 

7260  OUTPUT  709  ;"AC3VT3" 

7265  aAIT  100 

7270  ENTER  7o9  ;  0(1,2) 

7275  NEXT  I 

7285  CLEAR  (s  GCLEAR  £  DISP  "I  AM  FINDING  THE  ABSOLUTE  MA 
X  AND  FREQ  UPPER  AND  LOWER” 

7290  B 1  ,B4*C  @  H  ,Hl  ,H2  ,H  3  ,H4=  1 

7295  B2,E3*5U 

7300  FOR  1=2  TO  300 

7310  IF  D(I  ,2) <D(H,2)  THEN  7340 

7315  IF  D(I  ,2) =D(H ,2)  THEN  7325 

7325  Ao=D(I,2)  £F6=D(I,1)  &H*I 

7340  NEXT  I 

7345  A7=A6/SCR(2) 

7350  FOR  1*1  TC  H 

7355  IF  A7=D (1,2)  THEN  7405 

7360  IF  A7 <0 ( I  ,2 )  THEN  7365 

7365  IF  D ( I  ,2) <£ 1  THEN  7425 

7370  B  1=C  (1,2)  9  H  1  =  1  is  GOTO  7425 

7365  IF  C(I,2)>B2  THEN  7400 

7390  B2=D  (1,2)  ts  H2=I 

7400  GOTO  7425 

7405  El  ,B2=D  (1,2)  @  Hl,H2  =  I  9  F7=D(H2,j)  <H  GOTO  744u 
7425  NEXT  I 

7430  X*(A7-cl)/(B2-El) 

74  35  F7=X*  (D(H2,1).-C(H1,1))+D  (ril.l) 

7440  FOR  I =ri  TO  300 
7445  IF  A7=C (1,2)  THEN  7495 
7450  IF  A7 >D ( I  ,2 )  THEN  747  5 
7455  IF  C  ( I  ,2)  >B3  THEN  7470 
7460  E  3*C  ( I  ,2)  fc  H3*I 
7470  GOTO  7515 

7475  IF  C ( I  ,2) <*B4  THEN  74S0 

7460  B4*D(I  ,2)  t  H4  =  I 

74 9U  GCTO  7515 

7495  B3  ,£4*D (1,2)  H  H3,H4*I 

75u  5  L5*0  (H  3  ,  1) 

7510  GOTO  7530 
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7515  NEXT  I 

7  52G  X*  (A7-E4)  /  (E3-E4) 

7  525  L5«-  (X*  (D(H4  ,1).-D(H3  , 1)  )  ) +D  (H4  , 1) 

7530  G=F6/ ( L5-F7) 

7535  PRINT  "CENT EE  FREC  IS  ";F6 

7536  PRINT  "w  IS  ";C 

7538  DISP  "IF  YOU  APE  READY  TO  PROCEED ,  ENTER  'l'.  TC  RE 
RUN  FOR  BETTER  VALDES  ,  ENTER  '2''." 

7539  INPUT  Y 

7540  IF  Y*2  THEN  7235 
7545  RETURN 

blOO  REM  **CA1A  LIST* * 

8101  PRIN’IER  IS  2 

8102  FOR  1*99  TO  199 

8103  PRINT  USING  blC5  ;  D  (I  ,  1)  ,D  ( I  ,2)  ,D  ( I  ,  3 ) 

8105  IMAGE  DDCCCD  .C  ,2X  ,0  .DDDE  ,2X,D.DCCE 

8106  NEXT  I 

8107  PRIN1  USING  8  ICS 

8108  IMAGE  3/ 

8110  RETURN 

8200  REM  **PLGT  RE** 

8201  GC LEAR  @  CLEAR 

82G2  He*D  ( 1,1)  15*  (D(  3u0 ,1)-D(  1,1)  ) 

o2G  3  H9=C  (30C,1)  +  .1*(C  (  300 , 1  ).-C  (  1 , 1 )  ) 

8205  PLGTTER  IS  7u5  «  PEN  1 

8209  SCALE  He  ,H 9  ( . 8 *£ )  , £  +  . 2 5*E 

8210  XAXIS  0  ,  (C(  300  ,  l)-D<l,  1)  )/10  ,D  ( 3  ,  J)  ,D  (  300 , 1) 

8211  YAXIS  C( 1,L)  ,E/5  . 8*E)  ,E+.2*E 

8212  FEN UP 

8215  FOR  1=1  TC  300 
8  216  PLOT  0(1,1)  ,0(1,2) 

8217  NEXT  I 

8218  PE  NLP  y  DEG  t  LCIR  C,SIN(90) 

0221  FOR  LL  =  C  (1,1)  TO  0(  300  ,1)  STEP  (D  (  30U  ,  1 ) -C  ( 1 ,  1 )  )  /  10 
ig  FENUP 

8^22  mC v£  LI , - ( . 1 6*E ) 

8  22  3  LABEL  IN.T(Ll) 

8224  NEXT  LI 

6*.^b  LCIm  0  PENUP 

8227  MOVE  C  ( 1 , 1)  +  .2*  ID  (  300 , 1 )  -D  ( 1  ,  1 )  )  ,- (  .  7 * L ) 

8228  LABEL  "FREQUENCY  (H2)" 
o  2  30  PENuP  e  PEN  1 

o 2  .2  LD IR  0 

82  34  FOR  LL=- ( . 8 *E )  TO  E+ . 2*E  STEP  £/ 5  y  FENUP 
82  3  5  MOVE  D  ( 1 , 1 )  - .  09  *  (0  (  30C,1).-D  (1,1))  ,L  1 
82  36  LABEL  IN.1(L1) 

8  2  37  NEXT  LI 

8240  LCIR  C,SIN(9U) 

8  241  MOVE  C(i, 1  ).-.l*(0  (30b, L) -0(1,1))  ,.1*E 
8243  LABEL  ES 
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8250  RETURN 

8200  REM  **PLOT  IN** 

8  201  GCLEAR  g  CLEAR 

8  20  3  H8*D<  1,1).-.  15*  (D<  300,1) -D(  1,1)) 

8  304  H9*D  (300, 1)  +  .1*(D  (300, L) -0(1,1)  ) 

8305  L*MAX  (ABS(E5l)  ,ABS(E3)) 

8206  PLOTTER  IS  70S  @  PEN  1 
8212  SCALE  H8  ,H9  ,-L- . 5*L  ,L+ . 2 5*L 

6313  XAXIS  0  ,  ( C  (  30U  ,  1 )  -D  ( 1 , 1)  )  /  10  ,C  ( 1  , 1)  ,D(. 200,1) 

6214  YAXIS  D( 1,1)  ,L/5  ,-L-.2*L ,L+.2*L 

8315  PENOP 

6318  FOR  1*1  TO  20C 

8225  PLOT  C ( 1 , 1)  ,C(I,3) 

8330  NEXT  I 

8  2  21  PENUP  @  DEG  g  LCIR  0,SIN(90) 

6332  FOR  L1=C  (1,1)  TO  D(.2GG,1)  STEP  (D  (  30  C  ,  1 ) -D  ( 1 , 1)  )  /  10 
<§  PENUP 

6225  MOVE  L1,-L-.2*L 
6336  LABEL  IN.T(L1) 

8237  NEXT  LI 

6340  LOIR  0  6  PENUP 

6  341  MOVE  D(l,l)  +  .2*  (D(3CG,l)-0  (1,1)  )  ,-L-.4*L 
6342  LABEL  "FREQUENCY  (HZ)” 

6345  PENUP  g  FEN  H  LDIR  0 

6352  FOR  LL=-L-.2*L  TC  L-h.2*L  STEP  L/ 5  &  PENUP 
6255  MOVE  D(l,l)-.G9*fD(.3C0,l)-O  (1,1))  ,L1 
6  256  LABEL  IN-KL1) 
d  3  57  NEXT  LI 

6260  PENUP  £  PEN  1  £  LCIR  C,SIN(S0) 

8361  MOVE  D(l,l)-.  1*  (D  (  30G  ,  1 ) -D  ( 1 , 1)  )  ,~(L*.5) 

8368  LABEL  GS 

638G  RETURN 

8400  REN  **C IRCLES* * 

6401  GCLEAR  £  CLEAR 
8  <>02  L  1*MA X ( AES  ( E  3 )  , AES  ( E 5 )  ) 

6403  L*MAX  (AES  ( E )  ,L1) 

8404  PLOTTER  IS  705  £  PEN  1 

8406  SCALE  1. 25*  (-L- .  1*L )  ,  1 . 25*  (L+  .  1*L )  ,-L- .  2*L  ,L-K  .  15 *L 

8407  XAXIS  0,L/5.-L-.l*L,L+.l*L 
64G6  YAXIS  G  ,L/5  ,-L-.  1*L  ,L+  .  1*L 
8409  PENUP  £  PEN  1 

3411  FOR  1*1  TO  200 

6412  PLOT  L  ( I  ,2)  ,C<I,2) 

6413  NEXT  I 

6414  PENUP 

8415  LCIR  G  £  PEN  1  £  PENUP 
6416  MOVE  -  (L* . 5)  ,-L-.l*L 
c4l6  LABEL  H$ 

642G  PENUP  £  PEN  1  @  CEG 
6425  LCIR  0  , SIN ( 9G) 
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8426  FOR  LL»-L-.1*L  TC  L+.1*L  STEP  L/5  6  PENUP 

8427  MOVE  Li,-L 

8428  LABEL  INT(Ll) 

842S  NEXT  LI 

8435  LOIR  0  @  PEN  1  @  PENUP 

8436  FOR  LL*-L~. 1*L  TO  L+.1*L  STEP  L/5  g  PENUP 

8437  MOVE  -  (.886*L)  ,L  1 

8438  LABEL  IN.T(L1) 

8435  NEXT  LI 

8445  LOIR  G,SIN(9U) 

8446  fcCVE  -  ( . 95*L)  ,-(.3*L) 

8448  LABEL  G$ 

6450  RETURN 
8  5  45  REN  **NCRM* * 

6550  N8  ,N9=0 
8551  FOR  1=1  TC  10 

8557  WAIT  500 

8558  OUTPUT  709  ;"AC3" 

8559  ENTER  709  ;  N(I) 

8560  N8*N8+N(L) 

6565  OUTPUT  709  ;"AC4" 

8566  ENTER  705  ;  K(I) 

6  567  N5=NS+N(I) 

6568  NEXT  I 
6565  N8=N8/(I-1) 

6570  N9*-(N9/(I~1) ) 

8572  PRINT  USING  8573 
6  57  3  INAGE  3/ 

6575  RETURN 

8600  REN  ** VALUES** 

6601  CLEAR  BEEP 
8605  IF  L2=l  AND  M»2  THEN  66 1C 

6606  IF  L2=l  AND  N= 1  THEN  6611 

66  07  IF  L2=2  AND  N.  =  1  THEN  6b  12 

6608  IF  L2=2  AND  fc-2  THEN  6613 

8610  F6=F1  @  J=wl  @  GOTO  6615 

8611  F6=F0  s  w*w0  (?  GOTO  6615 

8612  F6=F6  £  W=W6  t§  GOTO  6615 

8613  F6*F9  e  w*w9 

ob  15  CISP  "PLEASE  SELECT  THE  HALF WIDTH  OF  INTEREST  TO  C 
LLECT  DATA.  (ENTER  HALF  WIDTH ) " 

ob  16  INPUT  Nts  CLEAR  t  CISP  "I  AN  CCLLECTING  THE  DATA." 
8625  FCR  1=1  TC  300 
8629  D  (I  ,1)=F6-N+I*N/150 
66  30  A5=" FR"  i  B 8  = " K 2 " 

6632  D$*v'AL$  ( C (  I  ,1)  ) 

6633  CLTPUT  717  ;A$,D$,E$ 

6634  WAIT  250 

6635  OUTPUT  709  ;"AC3vT3" 

6636  ENTER  70S  ;  K6 
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86  37  0(1  ,2) *H8*S8/N8 
8638  OUTPUT  709  ;"AC4VT3" 

8635  ENTER  70S  ?  H9 

8640  D  (I  ,3)*H9*S8/N9 

8641  NEXT  I 

8682  CLEAR  @  DISP  "I  AM  FINDING  MIN/MAX  VALUES  &  ASSOC. 
FREQS" 

6663  GCSUE  8700 

8664  IF  L2*2  THEN  8687 

86  65  PRINT  "GMAX*"  ;E  ,  FGMX=";E1 ,  , " EGMX*"  ;£ 2  ,  ,"EMX«"  ;£3  , 
,"FBMX*"  ; E 4  ,  ,"BNI*";E5  ,  ,“  FBMI=";£6 
8686  GCTC  8697 

66  67  PRINT  "RMAX  =  "  ;E  ,  FRMX=  ” ; E 1  ,  ,"XRMX=”  ;E2  ,  ,"  XMX*"  ;E3  , 
,“FXMX*";E4,  , " XM I * " ; £  5  ,  ,"FXMI  =  ";E  6 

8697  DISP  "IF  ALL  IS  nELL ,  ENTER  ''l'\  ENTER  '2'  1C  RETAK 
E  DATA." 

6698  INPUT  Y 

6699  RETURN 

d700  REM  **MAXMIN** 

6703  H  ,H 1 ,H2=  50 

6705  E*D  (50  ,2  )  £  E 1 ,  E4  ,£6=D  ( 5  0  , 1)  @  £  2  ,E  3  ,E  5*L  1 5  C  ,  3 ) 

8706  FOR  1*51  TC  250 

6710  IF  D  (I  ,2)  <D  <H,2)  THEN  6725 

8711  IF  C(I,2)*D(H,2)  THEN  6720 

8720  £=D  (1,2)  £  E  1=D  (1,1)  «  E2=C(I,3)  £  H  =  I 

8725  IF  D ( I , 3 ) <D (H 1 , 3)  THEN  6740 

6  726  IF  D  (I  ,3)=D (Hi  ,3)  THEN  8  735 

6735  E3=D  (1,3)  &  E4=D(I,1)  €  H  1*1 

8740  IF  D  (I  ,3) >*L(H2, 3)  THEN  8751 

3745  H2*I  @  E  5=D  (1,3)  «>  E6.=C(I,1) 

8751  NEXT  I 

8760  J1=D(H,2)  £  02*(D  (H+ 1 , 2) -D  (H-l  ,2)  )/2  «  J3*(D(H  +  1,2) 
+D  (H  - 1 , 2 ).-2* J L)/2 
3762  F6*-(J2/(2*J3)  ) 

8765  A=J  1+02+Fb+J  3*F6  “2 
8  770  Fb=D  (H  ,1)+F6*  (D  ( ri  ,  1 )  -0  (K-  1  ,  D  ) 

8775  E 1*F 6  €  E=A 
o 766  a*AES(£3-E5) 
a  7&  5  IF  M»2  THEN  675  3 
6790  IF  L2* 1  THEN  6792 
6  791  C 3=C  @  GCTC  o7£6 
3752  D 1  =  C  ts  GOTO  6  7  56 
8/93  IF  L2= 1  THEN  8795 

8794  D4*L  «l  GOTC  8756 

8795  C2*C 

6796  PRINT  "Dl="  ;D1,  ,"D2»";C2,  ,"C3  =  M;D3,  ,"D^  =  ";D4 

8799  RETURN 

9000  REM  **CALC~A** 

6001  PRINTER  IS  701,76 
50o2  PRINT  «§  PRINT 
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9003  PRIN1  "VALUES  MEASURED  IN  AIR  FROM  ";SS;"  DATA" 

9004  PRIN’I 

SU11  PRINT  "BLOCKED  CAPACITANCE*  ";A£S(CQ);"  FARADS  AT  " 
;F ; "HZ" 

9012  PRINT  "ELECTRICAL  QUALITY  FACTOR  *  "  ;Q0 

9013  PRINT  "MECHANICAL  QUALITY  FACTOR  *  ";w6 

9014  PRINT  "ELECTRICAL  RESONANCE  *";F0;“  HE" 

SC15  PRINT  "MECHANICAL  RESONANCE  *  ";F8;“  HZ" 

9016  PRINT  "DYNAMIC  ELECTRCMECHAN ICAL  COUPLING  COEFFICIE 
NT  *  " ?  * 

9017  PRINT  "STATIC  COUPLING  COEFFICIENT*  ";K1 
9021  PRINT  "BLOCKED  RESISTANCE*" ; RG; “OHMS” 

9031  PRINT  USING  9032 

9032  IMAGE  2/ 

9050  RETURN 

9100  REM  **CALC-vv* * 

9101  PRINTER  IS  701,76 

9102  PRINT  @  PRINT 

9103  PRINT  "VALUES  fCR  MEASUREMENTS  IN  WATER  FRCm  "  ;SS; 

"  DATA" 

9104  PRINT 

Sill  PRINT  "BLOCKED  CAPACITANCE  *  ";A£S(CC);"  FARADS  AT 
" ;  F  ;  "HZ" 

9112  PRINT  "ELECTRICAL  QUALITY  FACTOR  *  " ; w 1 

9113  PRINT  "MECHANICAL  QUALITY  FACTOR  *  "  ;Q$ 

9114  PRINT  "ELECTRICAL  RESONANCE  *  ";F1;"  HZ" 

9115  PRINT  "MECHANICAL  RESONANCE  *  ";F9;“  HZ" 

9116  PRINT  "DYNAMIC  ELECTRCMECHAN ICAL  COUPLING  COEFFICIE 
NT  =  " ; K 

9119  PRINT  "FREQUENCY  CF  OPTIMUM  EFFICIENCY  *  ";a0;"  HZ” 

9120  PRINT  "MECHANICAL  POWER  UTILIZATION  FACTCP" ; AES (P6 ) 

9121  PRINT  "  STATIC  COUPLING  COEFFICIENT*" ; K 1 
*122  PRINT  "BLOCKED  RESISTANCE  ="  ;  R9  ,*  "OHMS  " 

9132  PRINT  USING  9133 

913  3  IMAGE  3/ 

9150  RETURN 

920U  REM  **CALC-A/M** 

9201  PRINTER  IS  701,76 

9202  PRINT  «  PRINT 

9203  PRINT  "VALUES  MEASURED  IN  AIR  FROM  ";S9;"  DATA" 

9204  PRINT 

9211  PRINT  "ELCCRED  INDUCTANCE*  ";AoS(LC);"  HENRIES  AT  " 
;F6;"hZ" 

9212  PRINT  "ELECTRICAL  QUALITY  FACTOR  =  " ; w0 

9213  PRINT  "MECHANICAL  wUALITY  FACTOR  *  ";w& 

9214  PRINT  "ELECTRICAL  RESONANCE  *";F0;"  HZ" 

9215  PRINT  "MECHANICAL  RESONANCE  *  ”;Fb;"  HZ" 

9216  PRINT  "DYNAMIC  ELECTRCMECHAN ICAL  COUPLING  COEFFICIE 
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NT  *  ";K 

$217  PRINT  "BLOCKED  RESISTANCE*  " ; R 
$2  31  PRIM  USING  92  32 
$232  INAGE  2/ 

$250  RETURN 

$3u0  REM  **CALC-K/K** 

$301  PRINTER  IS  701,76 

9302  PRINT  £  PRINT 

9303  PRINT  "VALUES  FOR  MEASUREMENTS  IN  *AT£K  FROM  ";S$; 
"  DATA" 

9304  PRINT 

S3 11  PRINT  "BLOCKED  INDUCTANCE  *  ";A£S(L0);"  HENRIES  AT 
";F6?"H2" 

$312  PRINT  "ELECTRICAL  uUALITY  FACTOR  *  ";w3 
$313  PRINT  "MECHANICAL  UUALI'IY  FACTOR  *  "?wS 
9314  PRINT  "ELECTRICAL  RESONANCE  =  “  ;F1;“  HZ'* 

$315  PRINT  "MECHANICAL  RESONANCE  =  " ;F9;“  HZ" 

$316  PRINT  "DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 
NT  *  " ; K 

$316  PRINT  " POTENTIAL  EFFICIENCY  *  ";P5 

931$  PRINT  "FRECOtNCY  OF  OPTIMUM  EFFICIENCY  *  ";niO;H  HZ" 

9320  PRINT  "MECHANICAL  PC*ER  UTILIZATION  FACTOR" ;A£S (P6 ) 

9322  PRINT  USING  9323 
9  32  3  IMAGE  3/ 

$350  RETURN 

$$$&  CISP  "THE  END" 

$$$$  END 
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